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This  final  report  summarizes  the  research  performed  by  S-CUBED  under  AFOSR 
Contract  F49620-84-C-0029  during  the  period  from  March  1,  1984,  to  July  31,  1988. 

•  Partial  support  for  portions  of  the  research  was  provided  by  the  Defense  Nuclear 

Agency  under  Contract  DNA001-84-C-0127.  The  Co-Principal  Investigators  for  the 
project  were  Dr.  G.  A.  Hegemier  and  Dr.  H.  E.  Read,  The  AFOSR  Contract 
Technical  Monitors  were,  initially.  Lt.  Col.  L.  D.  Hokanson  and,  later  on.  Dr.  Spencer 
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and  consulting  company,  called  ENDOCHRONICS,  Inc. 
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ADVANCED  CONSTITUTIVE  MODELING  OF 
PLAIN  AND  REINFORCED  CONCRETES 


Section  1 
INTRODUCTION 


1.1  BACKGROUND. 

During  the  past  several  years,  the  Air  Force  has  been  deeply  involved  in  an 
extensive  effort  to  develop,  and  assess  the  feasibility  and  relative  effectiveness  of. 
various  candidate  modes  for  basing  the  MX  strategic  weapons  system.  In  most  of 
these  basing  modes,  the  key  elements  are  large  reinforced  concrete  structures,  called 
protective  structures,  which  are  designed  to  protect  a  missile  from  the  shock  loads 
prescribed  by  the  design  attack  scenarios.  The  enormous  costs  involved  in 
constructing  the  large  number  of  such  structures  required  by  the  system  dictates  that 
their  design  be  not  only  safe,  but  cost-effective  as  well. 

In  the  event  that  the  enemy  threat  changes,  it  is  also  important  for  the  strategic 
system  designer  to  know  the  ultimate  hardness  of  the  concrete  protective  structures, 
so  that  the  survivability  of  the  system  with  regard  to  the  new  threat  can  be  readily 
assessed.  The  most  expeditious  and  economical  way  to  do  this  is  through  the  use  of 
validated  analytical  models  of  the  structure's  behavior.  There  is,  accordingly,  a  need  to 
have  reliable  analytical  models  that  can  predict  the  loading  environments  for  which 
complex  reinforced  concrete  structures  will  collapse,  or  incur  sufficient  damage  to 
render  them  functionally  inoperable. 

1.2  OBJECTIVE 

The  objective  of  the  overall  research  program  described  here  is  to  construct  an 
advanced,  nonlinear,  multi-axial,  non-phenomenological  constitutive  model  of  reinforced 
concrete  that  will  provide  simulation  accuracy  in  the  nonlinear  response  regime  that  is 
superior  to  existing  models.  The  term  'advanced  nonlinear  multi-axial'  implies  a 
model  that  provides  greater  accuracy  than  existing  models  in  the  inelastic,  nonlinear 
response  regime  and  for  arbitrary  paths  in  multi-axial  stress  or  strain  space.  The 
term  'non-phenomenological'  implies  a  model  that  is  capable  of  synthesizing  the  global 
properties  of  reinforced  concrete  from  knowledge  of  the  plain  concrete  and  steel 
properties,  the  concrete-steel  interface  properties,  and  the  geometry  of  the  steel 
reinforcement. 


A  non-phenomenological  reinforced  concrete  theory  requires,  as  input,  constitutive 
models  for  each  of  the  basic  constituents,  namely,  steel  rebar  and  plain  concrete. 
The  accuracy  of  the  resulting  reinforced  concrete  theory  depends  heavily  on  the 
accuracy  with  which  one  can  model  these  two  constituents.  The  constitutive 
properties  of  steel  are  well  known  and  can  be  adequately  represented  by  reasonable 
simple  elasto-plastic  models.  Plain  concrete,  on  the  other  hand,  is  one  of  the  most 
complex  structural  materials  in  current  use  and,  despite  numerous  efforts  in  recent 
years  to  mathematically  model  its  properties,  there  is  today  still  no  model  which  can 
adequately  describe  its  nonlinear  constitutive  properties  over  a  wide  range  of  behavior, 
including  damage  accumulation  and  post-cracking  response.  As  a  result,  the  major 
challenges  that  one  faces  in  attempting  to  develop  an  accurate  constitutive  theory  for 
reinforced  concrete  are  twofold:  namely.  (1)  development  of  a  mixture  theory  which 
accurately  accounts  for  the  interaction  between  the  rebar  and  plain  concrete;  and 
(2)  formulation  of  a  constitutive  theory  for  plain  concrete  which  adequately  describes 
its  nonlinear,  inelastic  properties,  including  damage  and  cracking,  for  arbitrary,  multi- 
axial  load  paths.  Accordingly,  two  major  areas  of  research  have  been  pursued  in 
parallel,  concurrent  efforts  under  the  present  program;  namely,  (1)  development  of  a 
mixture  theory  which  accounts  for  concrete- rebar  interaction  in  a  realistic  manner;  and 
(2)  development  of  a  plain  concrete  model. 

Detailed  documentation  of  the  research  conducted  during  the  first  three  years  of 
this  study  can  be  found  in  the  Annual  Reports  by  Hegemier,  Read,  Valanis  and 
Murakami  (1985,  1986,  1987).  Details  of  the  research  performed  during  the  past 
(fourth)  year  are  described  m  the  topical  reports  by  Hegemier,  Murakami  and  Sweet 
(1988)  and  by  Read  (1988).  The  major  accomplishments  made  during  the  course 
of  the  entire  program  are  summarized  in  this  report,  and  the  reader  is  referred  to 
either  the  Annual  Reports  or  to  the  various  published  papers  and  Topical  reports  for 
further  details  of  the  research  accomplishments. 
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A  list  of  the  publications  which  resulted  from  this  contract  is  given  in  Section  4. 

These  topical  reports  are  included  as  part  of  the  present  report  in  Appendices  A 
and  B. 


Sectran  2 

REINFORCED  CONCRETE 


2.1  THE  PROBLEM. 

The  nonlinear,  inelastic  behavior  of  reinforced  concrete  is  dominated  by  complex 
interactions  between  the  steel  rebar  and  the  concrete.  These  interactions  have  a 
major  effect  on  structural  characteristics  such  as  stiffness,  strength,  damping  and 
ductility.  As  a  result,  it  is  necessary  that  a  model  of  reinforced  concrete  reflect  these 
phenomena  in  an  accurate  and  realistic  manner.  Further,  in  an  effort  to  minimize  the 
number  and  types  of  tests  necessary  to  define  model  parameters,  it  is  desirable  that 
the  model  be  non-phenomenological.  i.e.,  that  the  global  properties  of  reinforced 
concrete  be  synthesized  from  the  properties  of  (1)  the  steel  and  concrete,  (2)  the 
steel-concrete  interface  physics  and  (3)  the  steel  layout. 

2.2  ACCOMPLISHMENTS. 

A  theoretical  formulation  of  the  type  outlined  above  constitutes  a  new,  advanced 
model  of  reinforced  concrete  which  has  predictive  capabilities  far  superior  to  existing 
models.  In  an  effort  to  bring  such  a  model  to  fruition,  the  following  tasks  were 
accomplished  as  part  of  this  program: 

1.  Formulation  of  a  methodology  for  constructing  an  advanced 
theoretical  model  of  reinforced  concrete  which  correctly  reflects 
steel-concrete  interaction; 

2.  Construction  of  a  first  generation  model  of  reinforced  concrete  using 
the  results  from  the  above  task; 

3.  Construction  of  numerical  algorithms  and  a  special-purpose  computer 
program  which  allows  one  to  exercise  the  developed  model  for  a 
limited  class  of  test  problems: 

4.  Simulation  of  special  problems  (case  studies)  using  the  special 
purpose  code  developed  under  the  above  task. 

5.  Conduct  of  a  detailed  literature  search  and  evaluation  of  available 
experimental  data  on  the  behavior  of  reinforced  concrete. 
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6.  Performance  of  model  validations  consisting  of  experimental- 
theoretical  comparisons  of  important  response  features, 

7.  Performance  of  a  parametric  study  in  an  effort  to  determine  the 
influence  of  basic  material  and  geometric  properties  on  damage 
accumulation  and  failure  conditions. 

Under  Tasks  1  and  2.  a  new,  advanced  model  of  reinforced  concrete  with  a 
dense  unidirectional  steel  layout  was  constructed  (Hegemier  and  Murakami,  1985; 
Murakami  and  Hegemier,  1986;  Hegemier.  Read,  Valanis  and  Murakami,  1986).  The 
construction  technique  was  based  upon  the  use  of  multivariable  asymptotic  expansions, 
a  variational  principle,  and  certain  smoothing  operations.  The  resulting  model  was 
cast  into  the  form  of  a  binary  mixture  which  resembles  an  overlay  of  two  continua: 
steel  and  concrete.  These  continua  interact  via  body  forced  which  are  functionals  of 
the  relative  global  displacements  of  the  continua.  The  theory  is  fully  nonlinear  and  it 
incorporates  the  following  basic  physical  phenomena;  rebar  yielding,  steel-concrete 
bond  degradation  and  slip,  dowel  action,  and  progressive  concrete  cracking  (which  is 
treated  explicitly).  The  model  furnishes  both  global  and  local  (to  a  certain  degree  of 
accuracy)  measures  of  deformation,  stress  and  damage.  Response  characteristics  such 
as  stiffness  degradation,  ductility,  hysteresis,  strain  hardening,  and  certain  failure 
models  evolve  naturally  as  the  deformation  proceeds. 

The  validations  performed  under  Task  6  reveal  that  the  theoretical  framework 
developed  leads  to  a  model  of  reinforced  concrete  which  is  capable  of  accurate 
predictions  of  complex  response  characteristics.  (Hageman,  Murakami  and  Hegemier, 
1986;  Hegemier.  Murakami  and  Hageman.  1984;  Hegemier  and  Murakami,  1986; 
Murakami  and  Hegemier,  1986.)  Specifically,  the  validation  tests  performed  to-date 
clearly  indicate  that  the  model  correctly  simulates  progressive  concrete  primary  cracking 
(Hegemier,  Murakami  and  Hageman.  1984).  steel-concrete  bond  degradation  and  slip 
(Hageman.  Murakami  and  Hegemier.  1986;  Hegemier,  Murakami  and  Hageman.  1984), 
steel-concrete  dov\ic!  action  (Murakami  and  Hegemier.  1986).  and  certain  failure  modes 
(Hegemier,  Read,  Valanis  and  Murakami.  1986). 

Under  Task  7,  an  extensive  parametric  study  (Hegemier,  Murakami,  and  Sweet, 
1988),  which  is  included  herein  as  Appendix  A  for  the  convenience  of  the  reader,  was 
performed  using  the  developed  advanced  reinforced  concrete  model.  The  purpose  of 
this  study  was  to  ascertain  the  influence  of  fundamental  material,  interface,  and 


geometric  properties  on  the  nonlinear  response  characteristics  of  reinforced  concrete. 
Of  primary  interest  in  this  study  was  damage  accumulation,  overall  ductility,  and 
failure  modes/conditions.  The  damage  mechanisms  that  were  examined  included 
progressive  concrete  cracking,  steel-concrete  bond  slip,  and  rebar  yielding.  Parameters 
that  were  varied  included  the  concrete  tensile  strength,  the  steel-concrete  bond 
strength,  and  rebar  yield  stress  and  hardening  characteristics,  the  steel  volume 
fraction,  and  the  total  steel  surface  area.  Combinations  of  these  parameters  were 
discovered  that  allowed  a  concise  graphical  description  of  their  affects  on  material 
response.  The  latter  should  be  of  considerable  value  to  the  practicing  engineer. 

Finally,  during  the  last  research  period  the  unidirectional  theory  noted  above  was 
extended  to  include  multidirectional  steel  layouts.  The  latter  was  selected  in  the  form 
of  an  orthogonal  net.  Details  concerning  the  derivation  and  form  of  the  resulting 
theory  are  contained  in  Appendix  A. 

2.3  FUTURE  STUDIES 

It  is  recommended  that  future  studies  be  focused  upon  two  main  areas.  These 
are:  (1)  advanced  validation  tests  in  the  form  of  response  simulations  of  elementary 
structural  elements  under  shear,  bending,  direct  compression,  and  combinations  thereof; 
and  (2)  the  development  and  application  of  smoothing  operations  for  the  crack  field. 
Item  (1)  is  self-evident.  Item  (2).  however,  deserves  comment. 

Currently  the  model  development  process  homogenizes  the  steel  and  concrete  via 
appropriate  smoothing  operations.  These  operations  lead,  as  was  noted  previously,  to 
a  model  which  resembles  an  overlay  of  two  continua:  steel  and  concrete.  The 
evolving  concrete  crack  field,  however,  is  not  homogenized,  i.e.,  cracks  must  presently 
be  followed  explicitly.  This  feature  of  the  model  renders  it  difficult  to  use  from  a 
computational  standpoint.  As  a  consequence,  some  form  of  crack  field  smoothing 
should  be  explored  in  an  effort  to  simplify  the  model  computationally. 
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Section  3 

PLAIN  CONCRETE 


3.1  THE  PROBLEM. 

The  ultlir.jte  success  of  the  mixture  theory  for  reinforced  concrete,  discussed 
above,  hing-a  very  strongly  on  the  ability  of  the  constitutive  model  of  plain  concrete, 
used  with  the  theory,  to  accurately  describe  the  behavior  of  plain  concrete  over  the 
vide  range  of  response  that  can  be  expected  in  practice.  Any  deficiencies  in  the  plain 
concrete  model  will  surely  be  reflected  —  and  possibly  enhanced  —  when  the  model  is 
used  in  conjunction  with  the  mixture  theory  to  describe  reinforced  concrete.  Clearly, 
the  accuracy  of  the  mixture  theory  can  never  be  greater  than  the  accuracy  of  the 
plain  concrete  model  used  with  it.  Because  of  this,  every  effort  should  be  made  to 
develop  an  accurate  model  of  plain  concrete. 

During  the  past  ten  years,  considerable  research  has  been  devoted  to 
understanding  and  modeling  the  constitutive  behavior  of  plain  concrete.  From  this  has 
come  a  variety  of  different  constitutive  models  for  plain  concrete  which  provide 
reasonably  accurate  descriptions  of  plain  concrete  behavior  for  stress  paths  which 
generally  do  not  differ  greatly  from  the  standard  paths  followed  in  the  usual 
laboratory  testing.  Largely,  due  to  the  lack  of  available  appropriate  data,  little  proof¬ 
testing  of  the  models  has  been  done  to  examine  their  predictive  capabilities  for 
complex  stress  or  strain  paths  that  are  expected  to  occur  in  practice.  Many  of  the 
models  are  limited  to  the  stress  range  below  failure,  since  they  contain  no  provision 
for  treating  cracking.  A  few  of  the  models  do,  however,  attempt  to  treat  cracking 
and  post-cracking  behavior  as  well.  Generally  speaking,  however,  the  problems  of 
developing  a  constitutive  model  of  plain  concrete  capable  of  realistically  describing  the 
full  spectrum  of  behavior  from  pre-cracking  to  cracking  and  finally  post-cracking 
response  still  remains  basically  unsolved.  The  goal  of  the  present  research  is  to 
address  this  need. 

3.2  ACCOMPLISHMENTS. 

A  new  constitutive  model  of  plain  concrete  was  developed  which  appears  to  have 
remarkable  capabilities  for  predicting  the  nonlinear  inelastic  behavior  of  concrete  for 
stress  states  below  failure  (Valanis  and  Read.  1985;  1986).  The  model  is  formulated 
on  the  basis  of  the  endochronic  theory  of  plasticity  and.  as  such,  does  not  require  a 
yield  surface  nor  the  specification  of  loading  or  unloading  criteria,  as  in  classical 
plasticity.  It  predicts  that  plastic  flow  will  occur  from  the  onset  of  loading,  a  feature 
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which  makes  the  model  attractive  arid  appropriate  for  describing  the  behavior  of 
concrete,  which  does  not  exhibit  a  well-defined  yield  point.  Basically,  the  model  is 
isotropic  (when  referred  to  its  initial  state),  rate-independent  and  satisfies  the  second 
law  of  thermodynamics  (Clausius-Duhem  inequality).  It  realistically  portrays  the  major 
features  of  nonlinear  inelastic  behavior  exhibited  by  plain  concrete,  including  shear- 
volumetric  coupling,  effect  of  hydrostatic  pressure  on  shear  response,  hardening, 
hysteresis  and  stress-path  dependence. 

The  foregoing  model  was  applied  to  an  extensive  set  of  laboratory  data 
generated  by  the  University  of  Colorado,  using  a  true  triaxial  device.  The  test 
programs  consisted  of  six  different  series  of  stress-controlled  tests,  each  of  which  was 
designed  to  explore  a  particular  facet  of  material  behavior.  Altogether,  the  response 
of  concrete  to  over  45  complex  stress  paths  was  investigated.  In  all  cases,  the  stress 
paths  were  such  that  no  significant  macrocracking  occurred  during  the  tests. 

The  model  was  fit  to  a  very  small  subset  of  the  data,  after  which  it  was  proof- 
tested  by  driving  it  around  over  20  complex  stress  paths  to  predict  the  corresponding 
deformation  histories.  None  of  the  data  from  the  complex  stress  path  tests  were 
used  in  fitting  the  model  parameters  and  no  optimization  techniques  were  employed. 

As  shown  by  Valanis  and  Read  (1985,1986),  the  proof-tests  were  remarkably 
successful  and  revealed  the  powerful  predictive  capability  of  this  new  model.  In  all 
cases,  the  model  captured  the  essential  features  of  the  concrete  behavior  and  exhibited 
excellent  agreement  with  the  data.  For  nonlinear,  inelastic  behavior  inside  the  failure 
surface,  i.e.,  where  significant  macrocracking  has  not  occurred,  the  endochronic 
concrete  model  appears  to  provide  the  most  accurate  description  of  concrete  behavior 
of  any  of  the  existing  concrete  models. 

To  approach  the  problem  of  extending  the  basic  endochronic  model  to  account 
for  cracking,  anisotropy  and  dilatancy,  the  constitutive  behavior  of  a  brittle  elastic 
solid  was  first  considered.  From  this  effort,  there  resulted  a  new  continuous  damage 
theory  for  brittle  solids  (Valanis.  1985),  which  possessed  a  number  of  desirable 
features  not  found  in  other  models.  In  this  model,  damage  (microcracking) 
accumulates  in  a  gradual  manner  from  the  onset  of  loading,  provided  that  at  least  one 
of  the  principal  strains  is  extensional.  The  model  is  initially  isotropic,  but  becomes 
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anisotropic  if  damage  develops  in  specific  directions.  Through  a  clever  formulation, 
general  anisotropy  is  handled  through  the  use  of  second  order  tensors,  and  thus  does 
not  require  the  use  of  fourth-order  tensors,  as  do  most  such  models.  It  appears  that 
the  model  can  describe  the  standard  modes  of  cracking,  as  well  as  the  so-called 
"splitting  mode"  under  uniaxial  compression,  a  mode  which  is  beyond  the  scope  of 
most  existing  fracture  models.  Details  of  the  model,  together  with  a  number  of 
applications,  are  given  by  Valanis  (1985). 

Recently,  the  above  approach  for  describing  damage  and  fracture  of  solids  was 
extended  by  Valanis  (1987,  1988)  to  plastic-damaging  solids,  with  the  result  that  a 
new  endochronic-damaging  model  is  now  available  which  possesses  all  of  the  desirable 
features  found  in  the  earlier  (no  damage)  model  but  can  now  treat  cracking, 
anisotropy  and  dilatancy.  The  key  concept  behind  the  theory  is  a  mapping  which 
transforms  the  current,  damaged  and  generally  anisotropic  state  of  a  material  into  an 
isotropic,  undamaged  state.  The  model  describes  both  stiffness  degradation  and  yield 
limit  degradation  due  to  developing  damage,  as  well  as  dilatancy,  and  contains,  as 
special  cases,  several  models  that  are  noteworthy,  including  the  classical  elastic¬ 
fracturing  model.  The  model  has  been  successfully  apfjljed  to  recent  data  on  the 
response  of  plain  concrete  to  simple  tension  (Read,  1988)  .  For  further  details  of 

the  model  and  its  application  to  plain  concrete,  see  Valanis  and  Read  (1989). 

In  view  of  the  Air  Force’s  interest  in  shock  loading  of  defense  structures,  and 
with  the  goal  of  eventually  introducing  strain-rate  dependence  into  the  endochronic 
concrete  model,  a  comprehensive  review  of  strain  rate  effects  in  plain  concrete  was 
conducted  (Read,  1985;  Hegemier  and  Read,  1985).  In  particular,  we  sought  an 
answer  to  the  following  question:  Does  the  existing  data  on  strain  rate  effects  in 
plain  concrete  realistically  reflect  the  true  rate  dependence  of  this  material,  or  are  the 
inferred  rate  effects  the  result  of  spurious  system  effects,  inhomogeneous  deformation 
or  poor  methods  of  data  interpretation?  In  view  of  the  importance  of  strategic 
structure  response  to  this  study,  the  review  was  focussed  on  the  strain  rate  range 
from  10'^  to  10^  sec'^,  where  almost  all  of  the  data  comes  from  drop  hammer 
devices.  From  the  review,  it  was  concluded  that  considerable  caution  should  be 
exercised  in  using  such  data  to  develop  rate-dependent  constitutive  models.  The 
responses  of  test  specimens  in  such  devices  are  complicated  by  a  number  of  factors, 
including  inhomogeneous  deformation,  non-ideal  boundary  conditions  and  complex 
stress-wave  fields.  For  further  details,  the  works  by  Read  (1985)  and  Hegemier  and 
Read  (1986)  should  be  consulted. 


See  Appendix  B. 
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Concrete,  as  well  as  rocks  and  dense  soils,  when  compressed  at  constant  axial 
strain  rate  under  conditions  of  either  uniaxial  stress  or  triaxial  compression,  exhibit  a 
phenomenon  called  "strain  softening".  Materials  which  exhibit  such  softening  are 
characterized  by  a  constitutive  response  in  the  axial  direction  in  which  the  stress  rises 
monotonically  with  strain  to  a  peak,  and  then  decreases  with  further  increases  in 
strain. 

In  the  past,  strain  softening  has  been  generally  viewed  as  a  true  continuum 
material  property  and  routinely  incorporated  into  constitutive  models.  As  a  result,  the 
literature  abounds  with  advanced,  complex  constitutive  models  for  materials  such  as 
concrete,  rock  and  soil  which  are  designed  to  simulate  strain  softening.  Recently, 
however,  both  numerical  and  analytical  difficulties  have  surfaced  concerning  the 
solution  of  certain  wave  propagation  problems  in  strain  softening  materials.  As 
examples,  it  is  noted  that  strain  softening  can  lead  to  mesh  dependence  from  the 
numerical  viewpoint  and  loss  of  hyperbolicity  from  the  analytical  standpoint. 

In  view  of  these  difficulties,  and  because  plain  concrete  is  a  strain  softening 
material  (at  least  for  pressures  below  the  brittle-ductile  transition),  an  extensive  study 
of  strain  softening,  including  experimental,  theoretical  and  numerical  issues,  was 
conducted.  The  results  from  this  study  are  given  in  the  works  by  Read  and 
Hegemier  (1984),  Hegemier  and  Read  (1984),  and  Hegemier  and  Read  (1985).  In 
essence,  it  was  found  that  without  exception  the  initiation  of  strain  softening 
corresponds  to  the  transition  of  the  test  specimen  from  a  continuum  to  a  structure 
and/or  to  significant  geometrical  changes  in  the  specimen's  minimum  cross-sectional 
area.  It  was  concluded  that  strain  softening,  as  inferred  in  the  usual  manner  from 
conventional  laboratory  tests,  is  not  a  material  property  and  therefore  should  not  be 
incorporated  into  the  usual  (local)  types  of  constitutive  equations.  This  conclusion 
has  been  reaffirmed  by  subsequent  investigators. 

During  the  present  reporting  period,  there  has  been  considerable  activity  directed 
toward  developing  and  applying  non-local  models  of  strain  softening.  Virtually  all  of 
this  activity  has  focussed  on  the  case  of  softening  under  tension,  and  a  measure  of 
success  has  been  achieved.  The  more  difficult  case  of  softening  in  compression  has 
received  little  attention  and  continues  to  be  an  unresolved  problem.  Further  research 
is  needed  to  develop  a  general  approach  which  can  deal  with  softening  under  both 
tensile  and  compressive  loading  conditions. 
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Finally,  the  types  of  experimental  data  normally  used  to  develop  constitutive 
models  for  plain  concrete  were  critically  assessed  as  part  of  this  study  (Hegemier  and 
Read,  1985).  Various  issues  pertaining  to  strain  softening,  strain  hardening,  failure 
states,  failure  modes  and  strain  rate  effects  were  considered.  Particular  attention  was 
given  to  the  effects  of  different  boundary  conditions  and  different  test  devices  on  the 
resulting  data.  Also,  the  manner  of  defining  failure  was  explored.  It  was  found  that 
failure  modes  are  strongly  influenced  by  a  number  of  factors,  perhaps  the  most 
important  being  the  test  boundary  conditions.  Considerable  scatter  in  data  was  found 
with  respect  to  failure  modes,  especially  at  low  confining  pressures. 


3.3  FUTURE  STUDIES. 


In  order  to  complete  the  research  undertaken  in  this  task,  and  to  accomplish  the 
ultimate  goal  of  developing  a  constitutive  model  of  plain  concrete  that  is  capable  of 
describing  the  response  of  this  material  to  shock  loading,  including  damage  and 
fracture,  the  following  tasks  are  recommended  in  a  future  effort: 


•  In  view  of  the  very  encouraging  progress  made  in  this  study  toward 
developing  a  theory  for  the  plasticity  and  fracture  of  concrete, 
further  effort  should  be  undertaken  to  explore  and  validate  the 
theory  under  more  general  loading  conditions,  such  as  uniaxial 
compression  and  triaxial  compression.  These  cases  involve  cracking 
patterns  that  are  considerably  more  complex  than  those  which  occur 
under  simple  tension,  and  thus  should  provide  insight  for  further 
model  development. 


•  The  success  of  the  new  endochronic  plastic-fracturing  model  in 
describing  real  physical  phenomena  depends  strongly  on  the  ability 
of  the  damage  evolution  equation  to  reflect  the  underlying 
micromechanical  damage  processes.  Additional  effort,  therefore, 
needs  to  be  devoted  to  setting  the  damage  evolution  equation  on  a 
firm  micromechanical  foundation. 

•  Further  investigate  the  still  unresolved  question  of  strain-rate 
dependence  of  concrete,  with  particular  attention  on  the  critical 
assessment  of  experimental  techniques  developed  by  the  defense 
community  during  the  past  several  years. 
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If  the  issue  of  rate-dependence  can  be  resolved,  introduce  strain-rate 
effects  into  the  endochronic  plastic  fracturing  model.  Validate 
model  against  appropriate  strain-rate  data  not  used  in  calibrating 
the  model. 

While  considerable  progress  has  been  made  recently  in  under¬ 
standing  "strain  softening",  there  are  still  many  outstanding 
questions.  Most  importantly,  there  is  still  no  generally  accepted 
method  for  taking  strain  softening  into  account  under  general 
loading  conditions.  This  area  is  in  need  of  immediate  attention. 
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ABSTRACT 


This  report  describes  an  effort  to  develop  continuum  mixture  models  with 
microstructure  for  reinforced  concrete  with  unidirectional  or  bidirectional  steel 
reinforcement  layouts.  In  addition  to  model  construction,  analytical  and  numerical 
examples  are  presented  in  an  effort  to  demonstrate  the  capability  of  the  theory  to 
simulate  stiffness  degradation,  ductility,  and  failure  conditions.  The  theory 
incorporates  progressive  concrete  cracking,  steel-concrete  bond  slip,  and  yielding  of  the 
steel. 
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SECTION  1 
INTRODUCTION 

•  This  report  deals  with  microstructural  mixture  models  for  reinforced  concrete 

systems  with  unidirectional  and  bidirectional  steel  layouts.  Sections  2-7  concern  the 
development  of  models  for  3D  reinforced  concrete  elements  and  2D  reinforced  concrete 
panels  —  all  with  unidirectional  rebar  layouts.  These  models  are  examined  from  the 
standpoint  of  simulation  capability  in  Sections  8-10.  Included  here  are  discussions  of 
0  stiffness  degradation,  ductility,  and  failure  conditions.  In  Section  11,  the  modeling 

effort  is  expanded  to  include  rebar  layouts  which  are  bidirectional. 

The  models  developed  incorporate  progressive  concrete  cracking,  steel-concrete 
bond  slip,  and  yielding  of  the  reinforcement.  The  models  are  nonphenomenological  in 
^  that  all  model  parameters  are  defined  once  the  properties  and  geometries  of  the 

material  components  are  defined. 
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Section  2 

FORMULATION  OF  THE  LINEAR  PROBLEM 


Consider  a  reinforced  concrete  element  with  a  uniaxial,  periodic  steel  layout.  At 
time  t  in  the  loading  process,  let  this  element  consist  of  a  collection  of  intact 
subelements.  Let  V  and  9V.  respectively,  denote  the  volume  and  boundary  surface  of 
a  typical  sub-element.  Thus,  if  crack  surfaces  penetrate  the  element,  they  are 
confined  to  9V. 

For  reference,  let  rectangular  Cartesian  material  coordinates  x^,  x^.  x.  be  selected 
with  Xj  in  the  axial  direction;  Figure  2.1.  For  notational  convenience.  (  a  =  1.2. 
will  denote  quantities  associated  with  material  a  with  a  =  1  denoting  steel  and 
0  =  2  denoting  concrete.  In  addition,  the  notations  (') .  =  9(’)/8x.  will  be  employed. 
The  usual  Cartesian  indicial  notation  will  be  adopted  in  which  Latin  indices  range 
from  1  to  3.  and  repeated  indices  imply  the  summation  convention  unless  otherwise 
stated.  Quantities  of  the  form  (')  and  (  )  denote  dimensional  and  nondimensional 
variables,  respectively. 

Under  the  premise  of  small  deformations  and  linear  elastic,  isotropic  component 
response,  the  basic  governing  relations  for  the  components  and  the  component 
interfaces  are: 


(a)  Equilibrium  equations 


(a) 
j  '  J 


(2.1) 


where  a.,  denotes  the  stress  tensor,  represents  a  body  force  per  unit  volume, 
and  is  the  volume  occupied  by  material  a; 

(b)  Constitutive  relations 


on 


v(«) 


(2.2) 


where  are  Lamp’s  constants,  ffj,  is  the  Kronecker  delta,  and  e..  is  the 

infinitesimal  strain  tensor; 
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Figure  2.1.  Reinforced  concrete  subelement. 
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(c)  Strain-displacement  relations 


(a)  .  1 

•j  -  2 


(2.3) 


(d)  Interface  continuity  conditions 


on  8 

'J  J  J'  J 


(2.4) 


where  8  denotes  the  interface  between  the  two  constituents  and  is  the  unit 
outward  normal  to  material  1  on  3,  which  satisfies  =  0;  * 

(e)  Appropriate  boundary  conditions  on  3V. 

Conditions  (a)  -  (e)  define  a  well-posed  boundary  value  problem.  However,  a 
direct  solution  of  this  problem  for  the  stress  and  deformation  fields  on  an  element 
domain  V  is  rendered  virtually  intractable  for  most  cases  by  the  many  interfaces  3 
and  the  three-dimensional  aspect  of  the  problem.  A  primary  objective  of  the  work  to 
follow  is  to  alleviate  this  difficulty  by  a  process  of  homogenization  or  smoothing  of 
the  original  heterogeneous  material  in  V. 
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Section  3 

HOMOGENIZATION  PROCEDURE 


3.1  SCALING. 

The  homogenization  process  begins  with  a  simple  scaling  of  both  dependent  and 
independent  variables.  This  step  is  facilitated  by  the  introduction  of  the  following 
nondimensional  quantities; 


X .  =  X .  /A  , 
I  r  ' 


I  I  ' 


ij  ij  •  (m)' 


(3.1) 


(X, =  CX,J)(“>/E(^), 


where  A  denotes  a  reference  macro  length  and  E,  .  is  a  reference  composite  elastic 
modulus.  Under  (3.1),  equations  (2.1)  -  (2.4)  becorrie 


.  ft'') 

J  I 

i 

1 

u 

1 

IJ 

21  i,j 

(2) 

u>  ’  , 

U 


(3.2) 

(3.3) 

(3.4) 

(3.5) 


where  and  9  is  again  the  component 


iiucr 


A-9 


SSS-DFR-89-10119 


3.2  MICROCOOROINATES. 


It  is  expected  that  the  stress  and  deformation  fields  will  vary  significantly  with 
respect  to  two  basic  length  scales:  (1)  a  "macro"  length  typical  of  the  body  size  or 
loading  condition;  and  (2)  a  "micro"  length  typical  of  a  "cell"  planar  dimensions  as 
depicted  in  Figure  2.1.  These  macro  and  micro  scales  will  be  associated  with  the 
variables  A  and  A,  respectively.  Further,  it  is  expected  that  these  scales  will  differ  by 
at  least  one  order  of  magnitude  in  most  cases.  This  suggests  the  use  of 
multivariable  asymptotic  techniques.  Hegemier,  et  al  (1979).  This  approach 
commences  by  introducing  new  independent  variables  according 


^  0  as  e  ■*  0 
1>(0  1  as  e  -  0  , 

where  7  =  2,3  and 
e  5  A/A 


(3.6a) 


(3.6b) 


For  the  present  analysis,  it  will  suffice  to  set 
m  =  e  ,  |t(£)  =  1  . 


(3.7) 


Thus,  ail  field  variables  f(x.^t)  are  now  functions  of  the  "microcoordinates"  x*  as  well 

as  the  "macrocoordinates"  x^  =  x^,  i.c., 

7  7 


♦(‘i)  = 


(3.8) 


Spatial  derivatives  of  a  function  f(x.)  then  take  the  form 


f  .  =  F  .  +  e’*F  .• 
.1  .1  .1 


(3.9) 


where  ()  h  8(  )/9x,  and  (  )  h  a(  )/dXj  with  (  )  ,.  =  =  0- 

notationa'I  convenience,  the  furictions  f  and  F  will  both  be  written  as  f  i 
following. 


For 
in  the 


A-10 


SSS-DFR-8910119 


1 

•  3.3  SYNTHESIZED  FIELD  EQUATIONS. 

j 

The  operations  (3.8)  and  (3.9),  when  applied  to  all  field  variables,  lead  to  the 
following  synthesized  governing  field  equations  in  nondimensional  form: 


• 

♦  f(“^  =  0 

(3.10) 

J'.J  e  J'.J  ' 

ij  ij  xk  ij 

(3.11) 

• 

,(»)  ,  1  {(„(.)  .  „(«))  ,  ILW..  „(«).]), 

IJ  2  U  i,j  j,iJ  el  i,j  j,i  ))• 

(3.12) 

„(1)  .  „(?)^(1)  =  T-  on  8 

1  I  '  Jl  J  J‘  J  • 

(3.13) 

3.4  LOCAL  PERIODICITY  CONDITION. 

An  important  premise,  called  the  local  periodicity  condition,  is  now  introduced. 
This  condition  consists  of  the  assumption  that  local  periodicity  in  the  microvariables 
X*,  7=  2,3,  may  be  invoked  for  all  field  variables.  This  premise  allows  one  to 
analyze  a  single  cell  in  an  effort  to  determine  the  distribution  of  any  field  variable 
with  respect  to  the  coordinates  x*.  A  typical  such  cell  is  illustrated  in  Figure  3.1. 
The  local  periodicity  premise  is  suggested  by  the  Floquet  Theory  associated  with 
linear  differential  equations  with  periodic  coefficients,  the  set  of  relations  (3.10)  - 
(3.12)  can  be  reformulated  as  such  a  set  of  differential  equations. 

The  local  periodicity  condition  implies  that  all  field  variables  f(x.  ;x*)  satisfy 

K  7 


f[x^;  x*,x*;€)  =  f[x^;  -x*rx|;€) 


(3.14) 


where  (x*,Sp  and  denote  reflected  points  on  the  boundary  dtl  of  a  typical 

cell  with  area  Q  =  A^^uA^^ 
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Now,  on  the  boundary  0Q,  it  is  evident  that  Thus, 


X*  -  periodicity  implies  that 


(x^;x*;e]t/p)d(an)  =  0 


(3.15) 


dQ 


for  any  field  variable  f.  The  property  (3.15)  will  be  useful  later. 

3.5  WEIGHTED  RESIDUAL  PROCEDURE. 

In  this  subsection,  a  weighted  residual  procedure  is  introduced.  This  procedure 
will  be  subsequently  used  to  eliminate  the  microcoordinates  x*(7  =  2.3)  from  all  field 
variables  through  an  averaging  operation,  and  to  establish  appropriate  field  equations 
for  the  resulting  averaged  fields. 

To  begin,  let  ^  denote  the  space  of  all  H '-functions  f(X|^;x*)  on  V  with  respect 
to  X|^(k  =  1-3)  and  on  Q  with  respect  to  x*  (7  =  2,3)  that  are  x*  -periodic  according 
to  subsection  3.4.  Any  function  u.  e  #^with  u.  =  0.  on  8V\  where  u.  is  the 
specified  displacement,  will  be  called  an  admissible  trial  displacement.  Any  function 
5u.  £  f  with  5Uj  =  0  on  wiM  be  called  a  weighting  function  or  are  admissible 
variation  (of  u.).  In  addition,  let  f  denote  the  space  of  all  L^-functions  h(Xj^;x*)  on  V 
with  respect  to  X|^(k  =  1-3)  and  on  0  with  respect  to  x*  (7  =  2,3)  that  'are  x*- 
periodic.  Any  function  ffj.  e  f  with  =  f.  on  ^V.|.,  where  f j  is  the  specifiJd 
traction,  will  be  called  an  admissible  trial  stress. '  Any  function  6a..  e  #  with  6a..  =  J 
on  6\/j  will  be  called  a  weighting  function  or  an  admissible  variation  (of  9..).  In  the 
above,  it  has  been  tacitly  assumed  for  convenience  that  the  boundary  data  on  6V  is 
of  the  non-mixed  type,  i,e„  that  u.  is  specified  on  while  \.  is  specified  on  6V.j. 
with  5V  =  5V  and  5V  n9V.j.  =  0  (generalization  to  the  mixed  case  does  not 

represent  a  sigmficant  problem']. 
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» 


Next,  consider  the  weighted  residual  R  defined  by 
2 


R  = 


1 


a=l 

tfn  f  f 


du 

.(a) 


*  f.l 

'  I-  J',J  e  J|,J  IJ 


(a) 


dA 


(a) 


dV 


I 


a=l 


)  -  ..j)(">dA(“> 


,(«) 


dV 


i  .Tr[„p)  -  „p))d,. 


dV  ♦ 


I 


a=l 


(3.16) 


a'-) 


dS 


> 

0=1  J 


^  tH^i  -^j|«'j]^“^«uf“)ds’ 


dV  ♦ 


I  i[-a|Pi/j2))6u(2)d(8n) 

rao 


dV 


where  dV  =  dx^dx^dx^;  dA^*^  =  dx*dx*;  ds*  is  a  differential  element  of  0  £ 
dS  is  a  differential  element  of  aV;  atf  denotes  the  outer  cell  boundary;  T*  £  o, 
on  0;  for  a  given  trial  displacement  u,^®^  e.|“*  is  computed  according  to  (3.izj;  and 
the  notation  e..  (•••)  denotes  the  inverse  or  (3.11)  which  is  to  be  computed  given  a 


•i. 

trial  stress  a...  i.e., 
>j 


,(?)(...) ,  i-t  ,(«) .  idfl  „(»)« 

'  ‘  e(«)  ij  eW  "kk  ®ij 


u 


(3.17a) 


where 


f(o)  =  i«^®^(3X  > 


(X  ♦  2fi) 


(a) 


(«) 


2(  X  ♦  /i)^®^ 


(3.17b) 


In  (3.16),  the  integrations  with  respect  to  the  macrocoordinates  are  carried  out  over 
the  entire  domain  V  while  those  with  respect  to  the  microcoordinates  x*  are 
performed  over  a  typical  cell. 

^  satisfied  for:  (1)  all  admissible  iJu!®^  which  are  arbitrary  over  V 

and  A  ,  and  on  aV.^.  and  A'®^;  (2)  all  admissible  which  are  arbitrary  over  V 
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and  and  on  6\f^  and  and  (3)  ail  admissible  tff.  which  are  arbitrary  over 

V  and  8.  then  the  first  five  integrals  must  vanish  independently  and  it  is  evident  that 
weak  solutions  of  the  local  field  equations  have  been  generated  that  satisfy  the 
equilibrium  equations  (3.10)  (first  integral),  the  constitutive  relations  (3.11)  (second 
integral),  the  interface  displacement  continuity  (3.13)  (third  integral),  the  specified 
traction  on  8V  (fourth  integral),  and  the  traction  continuity  on  8  (fifth  integral).  The 
kinematic  condition  (3.12)  is  satisfied  identically  since  e..  is  computed  from  u.  via 
(3.12).  The  last  (sixth)  integral  vanishes  if  the  stresses  are  required  to  be  in  a 
solution  set;  otherwise  the  vanishing  of  this  integral  imposes  the  x*-periodicity 
condition  on  the  transverse  stresses  p..  (i,j  /  1)  (note  that  =  o”^on  the  cell 
boundary  80). 


If  R  =  0  is  satisfied  for  all  admissible  6u.  in  a  subspace  of  f,  and  for  all 
admissible  6a..  in  a  subspace  of  fl.  then  the  field  equations  (3.10)-  (3.13)  are  satisfied 
in  a  weighter/^  residual  (approximate)  sense. 


From  (3.16)  with  R  =  0.  Gauss*  Theorem,  and  the  x*-periodlcity  condition,  one 


obtains 


2  r 


Osl  J 


,(•) 


e  T  "  ^^i^)  *  ^’i(‘'P^  ■  up^]]ds» 


dxjdxjdxs 


• 

> 

t 

a=l 

< 

f  ^^^Su^^^dxjdxi 
<  I  2  3 


dxjdx2dx3 
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I  I 

(a) 


dx*dx| 


dS 


where 


5u. 


i.J 


Su. 

J 


(3.18) 


(3.19) 


Equation  (3.18)  is  similar  to  Reissner's  mixed  variational  principle  (1984.1986)  in 
which  variations  of  the  stresses  are  considered  along  with  variations  of  the 
displacements.  The  terms  of  (3.17)  involving  and  6T*  are  constraint  conditions 

that  reflect  satisfaction  of  the  material  constitutive  relations  and  the  interface 
displacement  continuity  conditions,  respectively.  Consequently.  (3.18)  can  be  envisaged 
as  the  principle  of  virtual  work  for  the  synthesized  fields  with  constraint  conditions; 
here  and  dT*  play  the  role  of  Lagrangian  multipliers.  In  what  follows,  only 
variations  in  the  transverse  stresses  rr.j(i.j/l.l)  will  be  considered,  i.e..  da^^  H  0. 


3.6  TRIAL  FUNCTIONS. 


As  indicated  in  the  previous  subsection,  the  weighted  residual  procedure  (3.18) 
will  be  subsequently  employed  to  eliminate  the  x*-dependence  (i.e.,  the  dependence  on 
the  microcoordinates)  from  all  field  variables,  and  to  establish  field  equations  for  the 
x*-averaged  (i.e.,  cel  I- averaged)  fields.  To  accomplish  this  task,  however,  appropriate 


e  f  and  e  f  must  be  postulated  or  constructed.  In 


tPial  functions  u.’  _  ,  _  . ..  ,  _  _  _  _ 

i  ij  ' 

particular,  it  will  be  necessary  to  exhibit  an  explicit  x*-dependence  of  al[  displacements 
u!®^  and  the  transverse  stress  components  (i.j  ^  1,1). 


The  required  x*-dependence  of  the  displacements  and  the  transverse  stresses  can 
be  obtained  by  application  of  an  asymptotic  procedure  developed  by  Hegemier  (1974) 
and  Murakami,  Maewal  and  Hegemier  (1981)  (see  Murakami  and  Hegemier  (1986)  for 
application  to  a  similar  problem).  This  procedure  is  as  follows:  The  premise  that 
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the  ratio  of  micro-to-macrodimensions  is  small,  e  <<  1,  and  the  form  of  the 
synthesized  equations  (3.10)  -  (3.17).  suggest  an  expansion  of  the  dependent  variables 
in  the  asymptotic  series 


n=0 


(3.20) 


If  (3.20)  is  substituted  into  (3.10)  -  (3.13).  and  the  coefficients  of  different 
powers  of  £  are  equated  to  zero  in  the  usual  manner,  a  sequence  of  boundary  value 
problems  (called  microboundary  value  problems  or  MBVP’s)  are  defined  on  the  cell. 
In  general,  solution  of  this  sequence  is  difficult.  However,  for  those  cases  where  a 
circular  cylinder  approximation  of  the  cell  boundary  5Q  is  appropriate  (i.e.,  where  the 
cell  aspect  ratio  is  approximately  unity),  consideration  of  this  sequence  to  0(e) 
motivates  the  following  form  for  the  trial  fields: 

(a)  Displacement  trial  functions 

♦  e  ij(x,^)g^®^(r)sin0  ♦  o(e^],  (3.21a) 

where 

§2  =  ^3  .  (3-21b) 

.  (3.210 

Xj  =  r  cosd,  x|  =  r  sin5  (3.21d) 

and  where  n^®^  denote  volume  fractions,  i.e.,  n^®^  =  where  A  =  A^^^UA^^^  =  Q 

is  the  total  cell  area. 

(b)  Transverse  stress  trial  functions 
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(a) 

(0) 

°22 

T22 

‘^33 

(’‘k'’'r'}  = 

“^33 

‘^23. 

.‘^23. 

'a2 


‘22’  (’■k) 


cos29 

cos2d 

0 


*  *3?  '‘k) 


cos2d 
-cos2d 
s  i  n20 


•  ‘if  (■'k) 


s  i  n2ff 
sin2d 
0 


•  I  '■2W3‘“'('') 


3cosd 

COS0 

Isind 


♦  f '’3(«k)9‘°’« 


s  i  n0 
3sin0 

COS0 


^31 

(«) 

“^31 

(«) 

► 

sln2^ 

('k)  *  7" 

‘if  ("k) 

°\2. 

^12. 

.cos29. 

•  ‘if  (‘k) 


cos2d 

si  n0' 

.-sin20. 

cosB. 

.  o(^2) 


(3.22a) 


(3.22b) 
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SECTION  4 

DEVELOPMENT  OF  3D  MIXTURE  MODEL  FOR  R/C 

Following  substitution  of  the  trial  fields  (3.21)  and  (3.22)  into  the  variational 
principle  (3.18),  integration  over  x*.  and  appropriate  integrations  by  parts,  one  obtains 
the  Euler- Lagrange  equations  of  the  variational  principle  in  the  form: 

(a)  Mixture  Equilibrium 


P.  .  =0  (i  =  1-3) 

JI.J  '  I  I  ^ 


(4.1a) 


where 


U 


1 

A 


M) 


C4.1b) 


cA 


P-  5  rr  i  Tfds*  = 


1_ 

cA 


(p  ds* 

a 


(4.1c) 


•  In  (4.1a),  is  a  cell-averaged  stress  tensor  and  P.  represents  an  effective 

interaction  body  force  due  to  the  transfer  of  the  stress  vector  across  the  interface  8. 
For  the  circular  cylinders  approximation.  A  =  r. 


(b)  Mixture  Equilibrium  -  stress  moments 


^(la) 

A2») 

1 

ffe(2)  ^ 

=  0  ' 

<^22 

-  ^22 

'arO) 

[^22 

,(2a) 

1 

f  t(2) 

•  m’]  =  0  ■ 

• 

^33 

‘'33 

2nW' 

1  ^22 

(4.2a) 


(4.2b) 
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,(la) 

23 

,(la) 

31 

,(la) 

12 


^(2a)  .  _±_  ^(2)  ^ 


23  23 


=  0 


.(2a)  ^  J_  A2)  _  0 

^^31  (1)  ^31  " 

n ' 


.(2a)  .  _l_t(2)  _  0 

^12  (1)^12  ~  °  ' 
n 


(4.2c) 
(4. 2d) 
(4.2e) 


(c)  Constitutive  Relations  -  transverse  stress  averages 


rx  ] 

-  u(“>  . 

r2XM- 

i  xH 

°Tl 

[PJ 

33  ^2,2 

1x2  J 

a+1 


y(®)  4.  Irll  §  33'\ 

'^1,1  (a)  ^2 


{\_±2ji] 


id) 


- 

"33 


U2. 


.0+1 


,0+1 


_J:_  ^(M)  .  u(«)  ,  *  2  i 

^(a)  "23  -  ^2,3  ^3,2  ^  Sj 

^(a)  "31  -  ^3.1  ^1.3  ^(a)  h  • 

_J_^(aa)  _  ^j(a)  ^  ^(a)  ^  § 

~U^12  -  ^1,2  ^2,1  (a)  ^1  • 


(4.3c) 

(4.3d) 

(4.3e) 


In  the  derivation  of  (4.3),  use  was  made  of  the  equality 
(1.1).  In  addition,  the  parameter  X  above  is  defined  by 

X  5  4/»(X  ♦  fi)  . 

(d)  Constitutive  Relations  -  transverse  stress  moments 

(2) 


for  (i,j)  ^ 


(4.3f) 


*(2)  .  .  SJL±Jil}y\l  .  §  1 
^22  ■  (2)  r2  ' 


(4.4a) 
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• 

‘33' '  ■  ^  (^2  •  y  ■ 

(4.4b) 

• 

[1^  .  ij  .  .  31^  1^, 

n  n 

(4.4c) 

t(2)  _  .  i£l  § 

^2  -  J2)  n 

(4.4d) 

• 

t  (2)  .  § 

^31  -  J2)  ^1  ■ 

(4.4e) 

(e)  Constitutive  Relations  -  Interaction  terms 

• 

(up).u«] 

(4.5a) 

e 

^2  -  P(2) 

(4.5b) 

• 

3  "  ^(3)  g2 

(4.5c) 

where 

• 

(4.5d) 

0=1 

• 

.-1  --1  X  h^®)a*3«)(“^ 

(4.5e) 

and  where 

• 

>»• .  i .  •»>  •  -“I 

(4.5f) 
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(f)  Constitutive  Relations  -  Axial  stresses 

(g)  Boundary  Conditions 

The  appropriate  boundary  conditions  for  (4.1a)  are 


either  =  0  or  j/. 

I  I  Ji  J 


(4.7a) 


on  dV,  where 


Y(ap)  =  1 


Yidx*dx|  =  i 


(a) 


Y.dxjdxl  . 


(4.7b) 


(a) 


The  relations  (4.1a),  (4.2a-e),  (4.3a-e).  (4.4a-e),  (4.5a-c),  (4.6),  and  (4.7a) 
constitute  a  complete  continuum  mixture  model  for  the  case  of  fully  three-dimensional 
fields. 
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SECTION  5 

DEVELOPMENT  OF  2D  MIXTURE  MODEL  FOR  R/C  PANELS 

In  this  section,  the  3D  mixture  model  of  Section  4  is  specialized  for  the  case  of 
plane  stress.  The  resulting  equations  constitute  a  model  for  the  in-plane  deformation 
of  reinforced  concrete  panels.  The  (material)  reference  coordinate  system  associated 
with  such  a  panel  is  depicted  in  Figure  5.1. 

5.1  PLANE  STRESS  CONSTRAINTS. 


The  following  constraints  are  deemed  to  be  appropriate  for  the  planar 
deformation  of  a  reinforced  concrete  panel: 


^(aa)  =  ^(aa)  ^  (aa)  ^ 
^^33  ^32  ^^31 


0,3  =  0 


U 


(a)  = 


=  0 


(5.1a) 

(5.1b) 


Equations  (5.1a)  and  (4.2b-d)  furnish 


.(2)  .  t(2)  =  0 

*31  "  *32  ■  ^  ' 


*(2)  _  2t(2) 

*22  "  “^*33 


(5.2a) 

(5.2b) 


The  relations  (5.2),  when  combined  with  (4.4)  yield 


I,  =  ij  .  0 

^3  '  [x  .'s/l] 

2i 
22 


(2) 


(la)  _  (2a)  _  A2). A\) 

°22  °22  ■  *22 


(5.3a) 

(5.3b) 


(5.3c) 


5.2  CONSTITUTIVE  RELATIONS. 

Under  the  constraint  (5.1a),  and  upon  noting  that 
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Material  coordinate  system. 
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i/E 


\  *  2ii  =  j  _  ^2  '  X  ♦  =  j  _  ^2  '  ^  =  2(1  +  I/) 


(5.4) 


one  obtains  from  (4.3a).  (4.3e),  and  (4.6) 


1  -  j/“ 


’■I 


.0+1 


^2,2"  (a)  ^2  "  ^1,11  ' 


.(aa)  _  r _ L_1 

‘^12  ~  12(1  ♦  u)} 


a+1 


(fl) ,  ,  izirii  I 


'1,2  ''2,1 


(a)  -1 


(5.5a) 

(5.5b) 


f(aa)  _ 
^11  ~ 


Aci) 


11  - 


"i!i  ^ 


u(fl)  ,  L-11°-  i  ' 

S,2  ^(0)  ^2 


(5.5c) 


The  expressions  (5.5a-c)  can  be  written  entirely  in  terms  of  and  as  follows. 
First,  the  quantity  can  be  determined  in  terms  of  u|®^  ana  by  substitution  of 
(5.5a,c)  into  (5.3c)  and  use  of  (4.4a)  and  (5.3b);  the  result  of  this  operation  is 


i  .  .  rj3!£_l  ‘'^u(l)  ,  f  T>E  f>  (2)  .  tE  1 


(1) 


u(‘) 

2,2 


(2) 


tE  ^‘'u(2 

il  -  u^}  2, 


(2) 

2 


(5.7a) 


where 


-1  _ 


w  E  1 

0)  1 

r  E  1 

(2)  2E^2^ 

ll  -  A 

2 

ll  - 1/^} 

*  n(«„ra(l  .  „W)(3  - 

Next,  substitution  of  (5.5b)  and  (4.4d)  into  (4.2e)  yields 


(5.7b) 
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where 


-1  e(»  F(g)fl  . 

.^(1))  *2n(*)„(2)(i  .  ,(2)) 


(5.8b) 


Finally,  upon  substitution  of  (5.7a)  and  (5.8a)  into  (5.5a-c),  one  obtains  the  desired 
constitutive  relations;  these  can  be  expressed  in  matrix  form  as 

W  =  |C){e}  (B.9a) 


where 

W.  „(2),(2»)T, 

<•>  ^  ■ 

[C]  5  ,  k,!  S  1  to  6  with  Cj^j  *  C||^  and 

C|^l  s  0  for  k  »  5,6  and  I  s  1  to  4;  (5.9b) 

and  where 

Cji  =  (n(')  -  ,  C,J  =  ,„(l)„WoO)p(2)  .  „Cl)„(2)c^^, 


*•22  ° 


„(2)  .  ^^(2)\(2)]„(2)  _  ^  („(1)  .  ^  pCl)j„(l)o(l)_ 


C„  =  70(«dW,.»)  =  ,  Cj3  =  70(l)D(«r(2)  =  .-Wc,,  , 

Cj,  =  ("®  -  toW)d®^W  ,  C3,  -  . 

Css  = 


fl-7  ^1 

C  -  n(2)„(2) 

1 .  2m  „(2)' 

V  n  j 

'  ^66  "  ^ 

[  nW'  J 

c  -  ,  .a)„(2)  r  -  „(2).{2) 

^56  ■  7(s)^  ^  '  ^66  "  "  ^ 


1  - 


^(2)1 

'^(s)n(2) 


(5.9c) 
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In  the  above. 


:(a) 


r(fl) 


D(a)  =  ■  ■■  Aa)  _  _ 

■  '2'/*  -2(1., (a)) 


1-1/ 


(a)‘ 


(5.9d) 


In  addition  to  the  material  constitutive  equations  (5.9a).  the  interaction 
'constitutive  relations"  are.  from  (4.5a.b): 


^  =  '’(1)  j  .  P2  =  ;> 


(of)  .  of)  ] 


(2) 


(5 . 10a) 


where  P^^yP^2)  expressed  in  terms  of  E.r/  by 


\l) 


aal 


eCa)  *  ^(2) 


h^'^^fl  ♦  i/^^“^f3  -  4»/) 


(«), 


(a) 


aal 


4E 


(a) 


(5.10b) 


5.3  EQUILIBRIUM  EQUATIONS. 

The  appropriate  equilibrium  equations  for  the  2D  model  can  be  obtained  directly 
from  (4.1a)  together  with  the  plane  stress  condition  (5.1a);  this  yields 

n(a)^(aa)  ^  0  (5.11) 

P1,P  77 

where  /J,7  a  1,2  . 

5.4  BOUNDARY  CONDITIONS. 

The  theory  is  completed  by  specification  of  the  boundary  conditions  on  9V;  these 
are,  from  (4.7)  given  by 

either  =  0  or  (5.12) 

P  Pi  1  P 

where,  again.  /7,7  =  1,2. 
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5.5  PRINCIPLE  OF  VIRTUAL  WORK. 


For  purposes  of  numerical  computations,  it  is  instructive  to  note  the  form  of  the 
principle  of  virtual  work  for  the  plane  stress  problem  under  discussion.  This  form  is 


[{a}^(5<e}  >  {F}'^5{U}]dV  = 


{T}'(5{U}dS 


(5.13) 


where 


<U}.{u(»,Up),U»),uM}T  _ 

<F} .  {.  Pj .  .  p^ .  Pj .  „rafm}T , 

<T}  3 

and  where  {e},  {a}  are  given  by  (5.9b). 

5.6  SUMMARY  OF  BASIC  EQUATIONS  IN  MATERIAL  COORDINATES. 

The  basic  equations  that  govern  the  linear  planar  response  of  a  R/C  panel  with 
a  unidirectional  steel  layout  are  summarized  below.  Use  of  these  equations  require 
that  the  reference  coordinates  be  material  coordinates,  i.e.,  that  the  x^-axis  be  aligned 
with  the  steel  as  depicted  in  Figure  5.1. 

(a)  Equilibrium  Equations 


where  a  =  1,2;  6,’]  =  1,2. 

(b)  Material  Constitutive  Relations 


(5.14) 


(ap)  > 


11 


uW  .  ^ 
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f(ap)  _ 
12  ■ 


^1,2  ^2,1  ^  (a) 


»<f)= 


u(«)  .  Jz.  ,  .WnW 

^2,2  ^  (a)  ^1,1 


(c)  Interactive  constitutive  Relations  for 
'“5  = 

where  6  =  1,2:  no  sum  on  S. 

(d)  Interaction  Constitutive  Relations  for  sf 

'  7 

$1  .  10(2)  (Uj  j  ,  J(2)  .  J(l)  , 


(e)  Parameters 

CW  ,  e(«)/(i  .  „2](«)  _  ;,(“).  e(“)/2(1  .  ./>(*)  , 

/’m  =  2  ^  .  ;f2,  =  -  4l/)(“)/8^(“) 

a=l  0=1 


h(‘)  =  1 


-1 


4^(2) 

,.(2) 


.(2) 


1  ■»  n 


(1)1 


.(1) 


(5.15) 


(5.16) 


(5.17) 


(5.18) 
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(f)  Partial  Stresses,  body  forces 


(5.19) 


(5.20) 


where  a  =  1.2;  6.j  =  1,2. 

(g)  Alternate  Form  of  Material  Constitutive  Relations 

M  =  lC]{e} 

where  {a},  {e},  and  [C]  are  defined  by  (5.9). 

5.7  BASIC  EQUATIONS  IN  GENERAL  COORDINATES. 


For  some  applications,  it  may  prove  convenient  to  utilize  reference  coordinates 
which  are  rotated  with  respect  to  the  material  coordinates  as  shown  in  Figure  5.2. 
The  basic  equations  can  be  rewritten  when  referenced  to  such  coordinates  as  follows: 

Let  x^.Xj  denote  the  material  coordinates  and  let  $^,X2  be  defined  according  to 


Xj  =  x^  cos^  -  Xg  sin^ 

A 

Xg  =  Xj  sin^  ♦  X*  cos^ 


(5.21) 


*3  "  ’‘3 


Then,  it  is  easily  demonstrated  that  the  stress,  strain  components  in  the  new 
coordinate  system  are  given  by 

2  -  2 

=  cos  ^  ♦  sin  ^  S22  -  2  sin^  cos^  Sjg  » 


’ll 


2  A  2  ^  ^ 

$22  =  sin  ^  Sjj  cos  ^S22  ♦  2  sin^  cos^  s^2 

*12  ~  '  sin^^]sj2  ♦  sin^  cos^[sjj  -  S22]  » 


(5.22a) 
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E 
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where  represents  or  the  latter  are  given  by 

=  U(“)  ♦  uW  (6,,  =  1,2). 


(5.22b) 


In  addition,  it  can  be  shown  that 
~  -  t2  sin^  , 

(5.22c) 

where  t_  represents  S^.  or  P  . 

7  111  1 

Substitution  of  (5.22a)  -  (5.22c)  into  the  basic  relations  of  subsection  5.6,  and 
with  use  of  (5.21)  to  evaluate  derivatives  with  respect  to  x^.x^  in  terms  of  derivatives 
with  respect  to  x^,X2.  one  obtains 

(a)  Equilibrium  Equations 

.  ;(.p)  .  0  (5.23) 

where  a  =  1,2  and  5,7  =  1,2. 

(b)  Material  Constitutive  Relations 

{a}  =  [C]{e}  (5.24a) 

where 

[c]  =  m  tc]  ;  (5.24b) 
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2 

cos  0 

0 

0 

2, 
cos  f 

•  2, 

sin  f 

0 

0 

•  2. 
sin  f 

2s i n^cos^ 

0 

•  2. 
sin  f 

0 

2. 
cos  f 

0 

-2sin^cos^ 

[T]  = 

0 

-si n^cos^ 

•  2, 

SI n  f 

0 

0 

s i n^cos^ 

2. 
cos  f 

0  1 

0 

2  2 
cos  ^-sin  ^ 

0 

-si n^cos^ 

0 

s i n^cos^ 

0 

The  6x6  matrix  [C]  is  defined  by  (5.9b.c). 
(c)  Interaction  Terms 


=  3^(0®  -  0<'>)  . 


where 

^11  ~  ^  ’ 

“  2  2 
P22  -  ^(2)*^°®  ^  "  ^(1)®'"  ^  ' 

^12  =  -  ^21  =  ('’(1)  - 


0 

2sin^cos^ 

0 

-2si n^cos^ 

0 

2  ■  2 
cos  ^-sin 

(5.24c) 


(5.25a) 


(5.25b) 
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SECTION  6 

MODIFICATION  FOR  INTERFACE  SLIP 

The  2D  and  3D  models  of  the  previous  sections  are  subject  to  the  restriction  of 
a  perfect  steel-concrete  interface  bond.  This  restriction  is  removed  in  this  section. 

Virtually  all  rebar  in  use  today  in  the  United  States  is  deformed.  Consequently, 
the  geometry  of  the  actual  steel-concrete  interface  is  complex.  For  modeling  purposes, 
it  is  conventional  to  replace  this  surface  by  an  artificial  mean  surface.  In  the  present 
analysis,  this  is  the  surface  d  in  dimensional  form  and  5  in  nondimensional  form. 

Given  the  mean  surface  9.  the  actual  steel-concrete  interaction  is  approximated 
by  specifying  an  interface  constitutive  law  between  the  traction  T*  on  9  and  the 
relative  slip  [u.]  =  u|^^  -  u.|^^  across  9.  This  law  must  reflect  local*  damage  to  the 
concrete  in  the  form  of  crushing  and  microcracking. 

For  most  practical  problems,  it  is  sufficient  to  incorporate  axial  slip  only.  The 
constitutive  law  associated  with  such  slip  can  be  expressed  generally  in  the  form 

T*  =  F([u^])  (6.1) 

where  F  is  a  functional  of  the  entire  slip  history.  A  subset  of  (6.1)  which  has  been 
previously  used  by  a  number  of  researchers  consists  of  a  slip  initiation  criterion  and  a 
slip  rule:  these  are  similar  to  the  yield  criterion  and  flow  rule  of  incremental  plasticity 
theory.  Application  of  incremental  plasticity  theory  to  slip  phenomena  has  been 
suggested  by  Drucker  and  Prager  (1952),  and  employed  by  Seguchi  et  at.  (1974), 
Bazant  and  Gambarova  (1980),  and  Bazant  and  Tsubaki  (1980). 

With  little  loss  of  generality,  one  may  assume  that  the  slip  is  axisymmetric;  see 
Murakami  and  Hegemier  (1986).  For  axisymmetric  axial  slip,  and  an  incremental  slip 
law,  the  phenomena  of  "slip"  may  be  incorporated  into  the  previously  derived  2D  and 
3D  models  by  replacing  the  constitutive  relation  for  P^  in  all  equations  with  a  relation 
of  the  form 
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(6.2) 


where  (  )  =  9()/9t.  t  denotes  (nondimensional)  "bookkeeping"  time,  and  is  a 


tangent  modulus  which  may  depend  on  the  field  variables. 


As  an  example,  the  following  slip  law  has  been  frequently  employed  for  cases  in 
which  reversed  slip  does  not  occur: 


=  ^1)  l'’ll<P(cr)' 


(6.3a) 


(cr) 


In  the  above. 


'’(cr)  '  2^^  "(cr)/' 


(6.3b) 


where  a.  .  >  0  is  the  (idealized)  local  critical  bond  stress.  The  law  (6.3a) 
corresponds  to  a  local  rigid-plastic  bond  law  of  the  form 

[uj]  =0  if  lT*j|  <  , 

[Uj]  ^  0,  T*j  =  (r^^^jsgn[uj]  if  )T*j)  = 


(6.3c) 


where,  again,  [uj  =  u|^^  -  u|^^  represents  the  local  interface  slip  across  9  in  the 
axial  direction. 

It  should  be  recalled  at  this  point  that  the  quantitites  [u.]  and  [U.]  measure 
different  phenomena.  The  former  on  9  means  local  slip  across  *9.  The  latter  (see 
(3.21a))  is  a  measure  of  the  difference  between  the  average  steel  and  concrete 
displacements  over  the  cell.  Consequently.  [U.]  includes  both  elastic  deformation  and 
slip;  the  quantity  is  a  measure  of  this  elastic  deformation. 
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SECTION  7 


MODIFICATION  FOR  REBAR  YIELD 


When  rebar  yield  occurs,  the  constitutive  relations  (5.9a)  must  be  revised.  If 
incremental  plasticity  is  utilized  to  describe  the  rebar  response,  then  (5.9a)  must  be 
replaced  by  an  incremental  relation  of  the  form 


ih  =  [DHe}  .  (7.1) 

The  matrix  [Dj  represents  a  tangent  modulus  matrix  whose  coefficients  may  depend 
on  the  plastic  component  of  {e}. 
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SECTION  8 

BASIC  EQUATIONS  FOR  UNIAXIAL  DEFORMATION 

A  problem  of  considerable  importance  concerns  the  simulated  behavior  of  a  R/C 
panel  when  the  latter  is  subjected  to  uniaxial  deformation.  The  relevance  of  this 
problem  stems  from  the  ability  to  conduct  laboratory  experiments  on  R/C  panels 
under  uniaxial  deformation  and  hence  to  validate  the  developed  model  by  comparison 
of  simulated  and  measured  response. 

The  uniaxial  problem  is  depicted  in  Fig  8.1.  The  panel  deformation  is 
constrained  to  the  x^-direction.  The  rebar  layout  makes  an  angle  ^  with  this 
direction.  The  governing  equations  for  the  R/C  panel  are  obtained  from  (5.23). 
(5.24a)  with  C,  C  replaced  by  D.D:  and  (5.25a)  with  replaced  by  Kj.i-  The 
appropriate  constraints  are 

=0  .  a  =  1,2  .  (8.1) 


Under  (8.1).  the  relevant  governing  equations  reduce  to 
(a)  Equilibrium 


I 


a 


(2p) 

11,1 


I 


where  (  )  j  =  d(  )/dXj  ; 


(b) 


Interaction 


^  2  2 
^11 "  *^(1)  ^  ^ 


(8.2) 


(8.3) 
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(c)  Constitutive 


s 

a 


5(1) 

iri,i 


D 


t 

a 


(2p) 

11 


B  .  D  u(2) 

‘^12^1,1  ‘^22^1,1 


(8.4) 


The  equations  (8.2)  -  (8.4)  are  employed  in  detailed  studies  of  simulation 
response  characteristics  in  the  following  two  sections. 
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SECTION  9 

STIFFNESS  DEGRADATION.  DUCTILITY,  FAILURE; 

UNIAXIAL  DEFORMATION  IN  DIRECTION  OF  STEEL  =  0) 

In  this  section,  the  nonlinear  response  of  a  R/C  panel  is  examined  in  detail  for 
the  case  of  monotonic  uniaxial  deformation  and  for  ^  =  0.  The  latter  condition 
implies  that  the  steel  is  aligned  with  the  direction  of  deformation;  see  Fig  8,1.  Using 
an  elastic-brittle  fracture  model  for  the  concrete,  an  elastic  plastic  model  for  the  steel, 
and  an  elastic-perfectly  plastic  description  of  the  steel-concrete  interface,  response 
characteristics  are  examined  that  relate  to  global  stiffness  degradation  of  the  R/C 
panel  under  extension,  ductility  of  the  R/C  panel,  and  failure  conditions.  The 
associated  analysis  incorporates  progressive  concrete  cracking,  steel-concrete  slip,  and 
yielding  of  the  rebar.  Of  particular  interest  is  the  influence  of  bond  strength  and 
concrete  cracking  on  the  overall  panel  response. 

9.1  FORMULATION  OF  THE  PROBLEM. 


The  formulation  begins  by  considering  an  initially  unloaded  R/C  panel  with  two 
starting  cracks  located  at  the  panel  termini  as  depicted  in  Fig.  9.1.  The  initial  crack 
spacing  is  It.  The  panel  loading  condition  corresponds  to  the  following  boundary 
conditions: 


0  at  a  £. 


u 

o 


at  X 


1 


e. 


(9.1) 


where  u  >  0.  The  notation  (  )  above  and  in  what  follows  indicates  that  the 
quantities  are  dimensional. 


Due  to  the  problem  symmetry,  one  need  only  consider  half  the  panel  as 
indicated  in  Fig.  9.2.  Upon  shifting  the  x^-axis  as  shown,  the  appropriate  boundary 
conditions  are 
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=  -  u  ,  =0  at  X.  =  0 

1  o'  11  1 

=0  at  Xj  =  1  . 


(9.2) 


Given  the  boundary  conditions  (9.2).  the  governing  equations  (8.2)  -  (8.4)  must 
be  solved.  Since  the  tangent  moduli  and  (k,£.  =  1,2)  depend  on  the  field 
variables  when  bond  slip  and  rebar  yielding  occur,  the  solution  for  a  particular 

loading  history  u  (t)  must  be,  in  general,  obtained  numerically.  However,  when  u^(t) 
is  monotonically  increasing,  an  analytical  solution  can  be  constructed.  For  this 
purpose,  the  domain  x.e(0.£.)  must  be  partitioned  into  appropriate  subdomains.  The 
vast  majority  of  the  desired  information  concerning  response  characteristics  can  be 
obtained  by  consideration  of  three  such  domains  as  follows  (see  Fig.9.2): 

Region  I:  Plastic  debonded.  In  this  interval  the  steel  is  plastic  and  the 
steel-concrete  interface  is  debonded. 

Region  II;  Elastic  debonded.  In  this  interval  the  steel  is  elastic  but  the 
steel-concrete  interface  is  debonded. 


Region  III:  Elastic  bonded.  -  The  steel  is  elastic  and  the  interface  is 
bonded  in  this  interval. 


Let  us  now  assume  that  (1)  the  steel  is  elastoplastic  with  linear  workhardening. 
Fig.  9.3:  (2)  the  concrete  Is  elastic-brittle  fracture;  and  (3)  the  interface  bond  is  rigid- 
perfectly  plastic.  Fig,  9.4.  Further,  in  an  attempt  to  simplify  the  analysis  with  little 
loss  of  pertinent  information  for  ^  =  0,  it  will  be  assumed  that  6^^  2  ® 

Under  these  conditions,  the  basic  equations  for  each  region  are: 


Region  I  (0  <  <  x  ) 


11.1  '  (cr) 


P' 

=  0  .  -  ?, 


'11.1  '  (cr) 


S  0 


11 


^(2P)  ,  „(2)|(2)ij(2) 

ail  "  ^  1,1 


(9.3) 
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where  P,  ,,  E,  ..  and  are  constants: 
(cr)  (p) 

Region  II  (Xp  <  <  xj 

^  =  0  •  -  Pr  ^ 

11,1  (cr)  11,1  (cr) 


=  0 


^(IP)  ,  a(2P)  =  n(2)^(2)g(2) 

^11  "  ^1,1  '  ‘^ll  "  ^1,1 

where  E(l)  is  a  constant. 


Region  III  (x 

1  K 

V 

<  0) 

®11,1 

V)! 

0j2)  .  =  0  , 

a(2p)  . 
®11,1 

'’(1)1 

[Op^  -  up^]  *  0  , 

)  = 

.®EWn(2) 

/  A 

(9.4) 


(9.5) 


9.2  SOLUTION  FOR  FIELD  VARIABLES. 

The  solution  of  the  governing  differential  equations  (9.3)  -  (9.5)  is  easily 
obtained  for  each  region.  These  solutions  can  then  be  pieced  together  at  the  region 
endpoints  by  imposing  continuity  of  both  displacements  and  axial  stresses 
Finally,  the  solution  is  completed  by  satisfaction  of  the  boundary  conditions  (9.2). 
The  result  is  as  follows: 

Region  I 
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(9.6) 


(9.7) 


cosh  0L(1  -  x) 
n^^^cosh/?£.[l  - 
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^(2a) _ Qnl 


il  ■  =  „(2)g(l)  "(O)*  l"(r,)»(s)  -  g(l)  ■ 

(9.8) 


a ,  ^x. 


coshglfl  -  x) 


.(1)  .  ^Cm)  ^(o) 


'11  -  „(2)|(2)g(l) 


a  sx, 


> 

V, 

cosh5£.fl  -  x) 

mi 

n(^)E^^)cosh  ^£.^1  - 

_ 

„(2)e(2)|(1) 


4> 


•  M  ~ 

\rx),^is)  -  ^(1)  ^(o) 


cosh/?g.fl  -  x) 


n(^)E(2)cosh  pt[i 


In  the  above,  the  following  variables  have  been  used: 


C 


2n(l)4‘’> 

~  -  X  Z  - 

^  =  l>^irxr  ' 


"(0)  '"u  (“)-  l'  %(1)J(»  *  „(2)f(2)  • 


(9.9) 


11  — 

(«  ^  '^  /^»)  ■  E(p)  =  plastic  modulus 


= 


f  '  1^- 

00 

L®/J 

1 

[2  .  nW  .  2  „(!)) 

/* 

(9.10) 


Further,  the  parameter  is  obtained  by  solution  of  the  transcendental  relation 


“W  *  0>t)tanh[/K.[l  -  *  E<« 


.Elaaflsl 


(9.11) 
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Overall  Displacement. 

In  addition  to  stress  and  strain  fields  in  each  region,  one  can  solve  for  the 
associated  displacements  0  Of  particular  interest  is  the  displacement  u,  .  =  - 
uj^^(O)  which  is  given  by 


(ml  - 

(rx)^s)  -  gd)  ^(o) 


'(o) 


,(i)e(i) 


tanh 


Pe-  '  „(2)|(2)|(1) 


(9.12) 


frxl 


f~2 

rc 


~2 


2n(l)E(l)  ^  (s)  (P)J 


]-? 


^(P))  • 


rx 


2n 


(Dp 


^(P) 


®(o)’‘(p) 


where 


(o)  -  p(D  F 

^  ^(P) 


=  Rebar  initial  yield  stress. 


9.3  FRACTURE  CRITERION  AND  SEQUENCE. 


(9.13) 


The  developed  mixture  model  of  reinforced  concrete  can  simulate  the  cell- 
averaged  fields  f^^^^(Xj^)  for  each  material.  The  model  can  also  be  used  to  estimate 
the  local  fields  f^^Hx|^;xy£)  within  each  material.  Consequently,  once  an  appropriate 
fracture  criterion  has  been  supplied,  then  one  can  proceed  to  examine  a  variety  of 
possible  fracture  modes.  A  subset  of  such  modes  include  primary  (i.e..  through) 
cracks,  secondary  (non-through)  cracks,  termini  "cone”  cracks,  and  axial  splitting 
cracks.  Although  all  such  fracture  modes  are  of  interest,  the  present  study  will  focus 
on  the  influence  of  primary  cracks  on  the  overall  response  characteristics  of  reinforced 
concrete. 
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An  elementary,  yet  useful,  criterion  for  the  initiation  of  a  primary  crack  can  be 
postulated  as  follows: 


-(2a)  x(c) 
max  ' 


(9.14) 


where  fj*"^  measures  the  "tensile  strength"  of  the  concrete. 

Examination  of  the  stress  field  in  each  of  the  foregoing  regions  reveals 

-  (2a)  - 

that  the  maximum  average  axial  stress  in  the  concrete,  max  a  *  occurs  at  x 

1*  1 


=  2- 


Consequently  in  the  absence  of  imposed  statistical  variations  of  f|^^  with  location,  the 
first  new  crack  (i.e.,  the  first  new  crack  system)  will  consist  of  a  single  crack  located 
at  the  mid-point  between  the  two  initial  termini  cracks.  The  value  of  the  applied 
stress  when  this  first  primary  crack  initiates  is  obtained  from  solving  (9.8)  with 
=  f|^^  and  x  =  1  together  with  (9.11).  Using  this  value  of  a,  ,,  the 


corresponding  displacement  u^^j  is  obtained  from  (9.12). 


(0)' 

Immediately  after  fracture, 
the  specimen  suffers  a  stress  drop  at  the  same  value  of  displacement  u^^j.  The  new 
applied  stress  corresponding  to  this  drop  is  obtained  by  replacing  t  by  £./2  in 

(9.12)  and  (9.11).  . 

additional  slip  occurs  during  the  fracture  process. 


It  is  noted  that  this  leads  to  jumps  in  a,  .  and  x, 

(rx)  (s) 


i.e.. 


Following  the  formation  of  the  first  new  crack,  the  R/C  element  consists  of  two 
subelements,  each  of  which  is  geometrically  identical  to  the  original  element  with  one 
exception:  The  subelement  length  (or  crack  spacing)  is  now  ^  rather  than  22-  Each 
subelement  is  now  extended  monotonically  from  the  value  of  the  displacement 
corresponding  to  the  initiation  of  the  first  fracture.  The  value  of  the  applied  stress 
at  a  given  displacement  for  each  subelement  is  given  by  (9.12).  (9.11)  with  2. 
replaced  by  2./2.  The  value  of  corresponding  to  the  initiation  of  a  primary  crack 
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within  each  subelement  is  obtained  from  (9.8),  (9.11)  with  x  =  1,  £.  replaced  by  2./2 
and  The  new  applied  stress  corresponding  to  the  stress  drop 

subsequent  to  fracture  is  obtained  by  replacing  £.  by  £./4  in  (9.12)  and  (9.11). 

Following  the  formation  of  the  primary  cracks  in  the  above  subelements,  the  R/C 
element  consists  of  four  subelements,  each  of  which  is  geometrically  identical  to  the 
original  element  except  that  the  subelement  length  is  now  £./2  rather  than  £_  Each 
subelement  is  now  extended  monotonically  from  the  value  of  the  displacement 
corresponding  to  the  initiation  of  the  second  crack  system,  and  the  entire  process  is 
repeated. 


The  foregoing  algorithm  leads  to  a  crack  sequence  as  depicted  in  Fig.  9.5.  A 
tyoical^  corresponding  stress-strain  curve  is  depicted  in  Fig.  9,6;  here  a = 
n'  ^fl'jl^^^(O)  represents  the  overall  applied  stress  based  on  the  entire  specimen  cross- 
section  and  e,  .  =  |u(o)/(£./N)|  is  the  overall  strain  (2£.  =  original  specimen  length,  N 
=  number  of  subelements). 

9.4  STIFFNESS  DEGRADATION. 

With  the  aid  of  the  crack  evolution  algorithm  described  in  subsection  9.3, 
stiffness  degradation  during  monotonic  extension  of  a  R/C  specimen  was  examined  in 
detail  for  the  case  ^  =  0,  i.e.,  when  the  direction  of  loading  was  aligned  with  the 
rebar.  The  "stiffness"  discussed  refers  to  either  the  tangent  or  secant  modulus 
associated  with  the  overall  stress  la,  J  versus  the  overall  strain  e,  see  subsection 
9.3. 

9.4.1  Influence  of  Bond  Strength. 

Bond  strength  was  found  to  have  a  major  influence  on  stiffness  degradation. 
This  influence  can  be  observed  in  Fig.  9.7  which  exhibits  predicted  relative  tangent 
stiffness  versus  global  strain  for  various  values  of  steel-concrete  bond  strength  ranging 
from  200  psi  to  •».  In  this  example,  the  steel  volume  fraction  n^^  =  0.01 
(1  percent),  the  concrete  cover  r  'Jli  =  0.10,  the  concrete  tensile  strength  f^  •  = 
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Figure  9.6.  Typical  stress-strain 


Relative  Stiffness  vs.  Global  Strain 

-  (7  <»)]/C7/»)  =  0.33,  =  400.0  psi 


Influence  of  bond  strength  on  stiffness  degradation. 
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400  psi,  and  the  rebar  hardening  characteristics  are  given  by  =  0.33: 

the  latter  is  perhaps  typical  of  grade  60  rebar. 

The  two  extreme  cases  in  Fig.  9.7  correspond  to  a  bond  strength  of  “  and  a 
bond  strength  of  200  psi  (strengths  in  the  range  600  -  800  psi  are  considered  typical 
of  many  practical  situations).  In  the  former,  stiffness  degrades  by  cracking  only  i.e., 
in  the  absence  of  slip.  In  the  latter,  stiffness  degradation  is  due  to  slip  only. 
Intermediate  cases  involve  both  slip  and  cracking.  It  is  evident  that  the  rate  of 
stiffness  degradation  increases  rapidly  with  decreasing  bond  strength,  i.e.,  increasing 
slip.  It  is  also  evident  that  tangent  stiffness  computed  in  the  absence  of  slip  (perfect 
steel-concrete  Interface  bond)  represents  an  upper  bound. 

9.4.2  Influence  of  Steel  Volume  Fraction. 

The  influence  of  steel  volume  fraction  on  stiffness  degradation  is  seen  by  an 
examination  of  the  200  psi  bond  strength  example  in  Figure  9.7  where  stiffness 
degradation  is  influenced  by  slip  only.  The  results,  for  steel  volume  fractions  ranging 
from  2  to  0.1  percent,  are  given  in  Figure  9.8.  The  rate  of  stiffness  degradation  is 
shown  '■Q  be  greater  for  smaller  steel  volume  fractions. 

For  cases  involving  both  slip  and  cracking,  the  same  trends  are  to  be  found 
(i.e.,  stiffness  degradation  occurs  more  rapidly  for  smaller  steel  volume  fractions).  It 
should  be  noted  that  the  sharp  decrease  in  overall  tangent  stiffness  for  steel  volume 
fractions  of  0.25  and  0.1  percent  in  Figure  9.8  is  due  to  the  rebar  undergoing  plastic 
deformation. 

9.5  DUCTILITY  AND  FAILURE. 

The  developed  mixture  model  of  reinforced  concrete  provides  stress  and  strain 
fields  for  each  material  component.  Consequently,  if  the  rebar  failure  condition  is 
specified,  then  the  mixture  model  can  be  applied  to  study  the  influence  of  various 
parameters  such  as  bond  strength,  steel  volume,  concrete  tensile  strength,  etc.,  on  the 
overall  ductility  of  the  composite.  Such  a  parametric  study  was  conducted.  The 
results,  which  appear  to  be  entirely  new.  are  presented  below. 
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Figure  9.8.  Influence  of  steel  volume  fraction  on  stiffness  degradation. 
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9.5.1  Derinition  of  Ductility. 

For  purposes  of  this  study,  a  ductility  measure  shall  be  defined  according  to 

=  e^/ej®^  (9.15) 

where  e|®^  denotes  the  failure  strain  of  the  rebar  alone,  and  e^  is  the  overall  failure 
strain  of  the  reinforced  concrete  composite.  Consequently,  =  1  for  the  rebar  alone 
while  <  1  for  the  composite. 

It  is  noted  that,  for  most  practical  situations,  the  strain  field  in  the  steel  will  be 
highly  nonuniform.  An  example  is  depicted  in  Fig,  9.9.  Consequently,  Q)  can  be 
expected  to  be  considerably  less  than  unity  for  most  cases  of  interest. 

9.5.2  Influence  of  Bond  Strength. 

Bond  strength  exerts  a  major  influence  on  overall  specimen  ductility.  This 
influence  is  exemplified  in  Figure  9.10  for  a  range  of  steel  percentages.  The  concrete 
tensile  strength  in  this  example  is  sufficiently  high  to  preclude  concrete  cracking  prior 
to  failure  of  the  rebar.  As  can  be  observed,  the  ductility  reduces  sharply  from  1.0  as 
the  bond  strength  is  increased  from  zero.  For  bond  strengths  in  excess  of  600  psi. 
the  ductility  is  generally  less  than  0.1  for  a  range  of  practical  steel  percentages  —  in 
the  absence  of  cracking. 

9.5.3  Influence  of  Concrete  Tensile  Strength. 

Concrete  cracking  —  or  concrete  tensile  strength  --  also  plays  a  major  role  in 
overall  specimen  ductility.  Figures  9.11  -  9.14  exhibit  the  influence  of  concrete 
cracking  for  a  range  of  steel  percentages.  It  is  evident  that  the  effect  of  progressive 
cracking  is  to  dramatically  increase  the  overall  ductility.  Cases  where  cracking  occurs 
are  indicated  by  the  curves  which  show  sudden  changes  in  S'.  A  comparison  of 
cases  reveals  that  ductility  in  general  increases  as  the  concrete  tensile  strength 
decreases. 

9.5.4  Influence  of  Rebar  Strain  Hardening. 

The  rebar  in  the  current  study  was  modeled  as  piecewise  linear  under  monotonic 

extension  with  an  elastic  modulus  of  3  x  10^  psi  and  a  plastic  modulus,  E  ,  of 

P 


Rebar  Strain 

r(2)/i  =  0-10.  nO)  =  0.50  (%) 

-  cry(“)]/(7y<*)  =  0.33,  (T<^)  =  400.0  psi 
x'/l  =  0.18,  x,/i  =  0.65 
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Figure  9.9.  Typical  nonuniform  strain  field  in  the  steel. 
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Figure  9.10.  Influence  of  bond  strength  on  specimen  ductility. 


Relative  Ductility  vs.  Bond  Strength 

[(J,«  -  ayW]/ay<»)  =  0.33,  r(2)/i  =  0.10 
fl')  =  200.0  psi 


SSS-DFR-89-10119 


(s)'a/'a  ‘iC;ini=>na  SATiBpa 


A-59 


500  1000  1500 

Bond  Strength,  (psi) 

Figure  9.11.  200  psi  concrete  tensile  strength. 
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Figure  9.12.  400  psi  concrete  tensile  strength. 
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approximately  2-3  orders  of  magnitude  less.  Within  the  context  of  this  description, 

the  ductility  was  found  to  depend  only  weakly  on  the  strain  hardening  modulus  E  . 

This  weak  dependence  can  be  observed  from  Figures  9.15  -  9.18  wherein  the 

modulus  E  varies  from  E  =  3  x  10  psi  to  E  =  6  x  10  psi. 

P  P  P 

9.5.5  Ductility  versus  Dimensionless  Bond  Strength. 

In  the  process  of  conducting  this  study,  it  was  discovered  that  one  can  collapse 
all  ductility  versus  bond  strength  data  by  graphing  that  ductility  versus  a 
dimensionless  bond  strength  defined  according  to  o  /E^  where  a  and  E^  are  given 
in  (9.9).  The  resulting  master  curve  is  shown  in^l^ig.  9.19.  Tfils  curve  applies  for 
cases  involving  cracking  as  well  as  for  those  wherein  no  cracking  occurs.  When 
cracking  occurs.  22.  represents  the  crack  spacing.  If  slip  only  occurs,  22.  is  the  length 
of  the  specimen. 
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Figure  9.15.  Rebar  plastic  modulus  E  =  3  x  10  psi. 
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Figure  9.16.  Rebar  plastic  modulus 
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Figure  9.17.  Rebar  plastic  modulus 
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Figure  9.18.  Rebar  plastic  modulus  E  =  6  x  10  psi. 
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r'gure  9.19.  Collapsed  ductility  versus  bond  strength  curve. 
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SECTION  10 

EXTENSION  OBLIQUE  TO  STEEL 

In  this  section,  the  response  of  a  R/C  panel  is  examined  for  the  case  of 
monotonic  extension  in  a  direction  which  is  oblique  to  the  steel  layout,  Fig.  8.1.  The 
focus  of  the  presentation  is  on  the  influence  of  the  angle  ^  on  stiffness  degradation. 
In  the  analysis  to  follow,  the  steel-concrete  bond  is  assumed  to  be  perfect. 

10.1  BASIC  EQUATIONS  AND  SOLUTION. 

The  governing  equations  for  this  case  are 

‘’u^l  •  =  °  ■ 

“IM  ■  (^(1)"“^*  • 

.  5  .  D 


J.(2p)  u  Jj(l)  Q  g(2) 


Substitution  of  (10.2)  into  (10.1)  furnishes 


A 

d^Gp) 

fl2 

*  ^11  • 

[Gp).Gp)] 

dxj 

=°. 

.  d^'GJ') 

n  ^ 

*“22 

d^(« 

J 

[Gp)-Gp)] 

dxg 

■  ^1 

(10.3) 


where 


^11  '  ^(1)*=®®^^  *  ^(2)®'"^^ 


(10.4) 
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and  where  p 


(2) 


are  given  by  (5.18). 


The  solution  of  (10.3)  is  readily  obtained  in  the  form 


^  =  Aje  ♦  AjC  ♦  AgXj  +  A^  , 


‘'(2)  ^*1  -Xx. 

Uj  =  -  c  A^e  -  cAje  ♦  *  % 


where 


°11  *  ^22  *  ^^12 

L^J 

.  ^11^22  "  ^12 

(10.5) 


c  = 


“fl  *  “l2(- 

»ii  *  “J 

.^11  *  ^12] 

[622  ♦  ^12 

' 

hz  •  “12] 

(10.6) 


Equations  (10.5)  lead  to  the  concrete  stress: 


'11 


,  n(2);(2a)  ,  5 

.  D  G(2) 

_  n  Pjj  -  1^12  1,1 

"22^1,1 

Xx.  “I 

A 

-Xx 

=  *1«  -  lOjjJ 

♦  Age 

*  *3K2  •  ^22) 

Djj  .  aojj] 


The  following  boundary  conditions  are  now  specified: 


(10.7) 
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=0  at  =  0  . 


=  £,  e  at  Xj  =  £.  ,  (e  specif  ied) 


=  0  at  x^  =  «.  . 


Substitution  of  (10.5)  and  (10.7)  into  (10.8)  furnishes 


(10.8) 


A  I 


Di2  ♦  D22]  =  - 


A  .  2e  a  r; 
''3  ■  *  A 


[Dj2  “  c  D22]  cosh  X£.  . 


(10.9) 


A.  =  0  , 


where 


A  =  2^Dj2  022]*' oh  X£.  -  2X£.^Dj2  ”  c  D22]cosh  X£. 


(10.10) 


The  solution  of  (10.3)  subject  to  the  boundary  conditions  (10.8)  can  now  be 
written  in  the  form 

0(1)  _  2i£:  ^  822)5 'nh  Xx^  -  x[o^2  *  ^  ^22)  ^^^[’'l)] 


(10.11) 


op)  = 


"  "  (^12  *  ^22)®'"^  ^’‘1  "  ^[^12  ■  ^  ^22]^^®®^ 


The  displacements  (10.11)  furnish  the  stresses: 


„(i),(i>) ,  Ip  .  022]  Ki  - '  ‘>12)'=“'’  ’“1 

-  x(Djj  -  c  D2j)(0j,  .  DjJcMh  Xi]  . 
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# 


n 


(2):(2a)  _  2e^fr; 
ni  ■  A  Lri2 


-  X  [Dj2  -  c  D22) [Dj2  ♦ 


-  c  D22jcosh  Xxj 
D22]cosh  X£. 


»  Examination  of  (10.12b)  reveals  that  attains  a  maximum  at  x^  =  0;  the 

corresponding  value  of  is  given  by 

n(2);(2a)  ^  ^  *  ^22^  [^12  '  ^22) 

Similarly,  examination  of  (10.12a)  reveals  that  max  occurs  at  x^  =  £.  where 

=  ^)  =  *'=)(^iA2  ■  ^12]*=°**’ 

10.2  FRACTURE  CRITERION  AND  SEQUENCE. 

Consider  now  the  problem  of  predicting  the  evolution  of  the  primary  crack  field 
in  the  R/C  specimen  when  the  latter  is  subjected  to  monotonic  extension.  If 
condition  (9.14)  is  again  adopted  as  the  criterion  for  primary  crack  initiation  in  the 
concrete  then,  since  max  again  occurs  at  the  specimen  center  x^^  =  0  (see 

(10.13)),  the  fracture  sequence  as  well  as  the  resulting  overall  stress-strain  behavior 
can  be  computed  using  an  algorithm  similar  to  that  outlined  in  Section  9.3. 

10.3  STIFFNESS  DEGRADATION. 

Figure  10.1  shows  the  degradation  of  overall  secant  stiffness  with  increasing 
overall  specimen  strain  for  rebar  angles  ranging  from  ^  =  0*  (rebar  along  loading 
direction)  to  ^  =  60*.  Remarkably,  the  degradation  curves  fall  into  a  reasonably  thin 
band  for  the  strain  interval  shown.  However,  while  the  stiffness  "decay  length"  is 
similar  for  different  rebar  angles,  major  differences  occur  in  both  the  "initial"  stiffness 


* 


Secant  and  tangent  stiffness  are  equivalent  when  the  interface  bond  is  perfect  and 
the  rebar  remains  elastic. 
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Relative  Stiffness  vs.  Global  Strain 

n(‘)  =  1.00(%) 

_ =  0.10,  fl'>=  400.0  psi 


Figure  10.1.  Degradation  of  overall  secant  stiffness. 
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(zero  strain,  E  1  and  the  "final"  stiffness  (•  strain,  E.,,  ).  The  variations  in  these 
initial  and  final  stiffness  with  rebar  angle  are  depicted  in  Figures  10.2  and  10.3  for 
steel  volume  fractions  ranging  from  0.1  to  1.0  percent.  Variations  in  the  final 
stiffness  are  masked  in  Figure  10.1  by  the  fact  that  they  all  tend  to  be  approximately 
an  order  of  magnitude  smaller  than  the  initial  stiffness. 

When  reviewing  the  data  depicted  in  Figure  10.1,  it  should  be  recalled  that  the 
steel-concrete  interface  bond  has  been  assumed  to  be  perfect,  and  the  rebar  has  been 
taken  as  elastic.  As  can  be  observed  in  the  examples  cited  for  the  case  ^  =  0*. 
steel-concrete  interface  slip  and  rebar  plasticity  can  be  expected  to  alter  the  foregoing 
results  to  a  significant  degree  in  some  cases. 


Initial  Stiffness  vs.  Rebar  Angle 

r<2)/i  =  0.10 


^Tsd^l  O'J  ‘SSaUIIIlQ  TBtlTUT 


Figure  10.2.  "Initiar  (zero  strain)  stiffness. 


Final  Stiffness  vs.  Rebar  Angle 

r<2)/i  =  0.10 


strain)  stiffness 
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SECTION  11 


EXTENSION  OF  MIXTURE  MODEL  TO  ORTHOGONAL 
BI-DIRECTIONAL  REBAR  LAYOUTS 

The  mixture  models  discussed  in  the  preceding  sections  apply  for  uniaxial  rebar 
layouts.  Most  practical  cases, however,  concern  bi-directional  or  tri-directional  layouts. 

In  what  follows,  the  previous  mixture  theory  is  extended  to  include  the  case  of 
bi-directional  orthogonal  rebar  layouts.  Both  rebar  plastic  deformation  and  steel- 
concrete  interface  slip  are  included  in  the  analysis. 

11.1  A  3D  CELL  FOR  TWO  REBAR  SYSTEMS. 

The  R/C  material  to  be  considered  is  depicted  in  Figure  11.1a.  As  indicated, 
two  reinforcement  systems  are  now  involved.  These  are  designated  as  Rebar 
System  1  and  Rebar  System  IL  Both  systems  are  presumed  to  be  periodic.  This 
premise  leads  to  the  3D  cell  illustrated  in  Figure  11.1.  This  cell  is  adopted  as  the 
basic  building  block  in  the  following  development. 

11.2  SCALING  AND  MICROCOORDINATES. 


With  reference  to  Figure  11.1,  the  quantity 


l/3/,l/2 


(11.1) 


is  selected  to  represent  the  cell  "size"  and  the  parameter  e  =  A/A  is  adopted  as  the 

micro-to-macro  dimension  ratio  where  the  meaning  of  A  is  as  was  indicated  in 

Section  3.1,  i.e.,  A  is  a  reference  macro  length.  Next,  the  space  variables  x  and  x* 
are  defined  according  to 

X  =  X  /A  ,  X*  =  X  /A  .  (11-2) 

As  in  the  previous  analysis  for  the  unidirectional  steel  layout,  the  components  x.  and 

X*.  represent  "macro"  and  "micro"  coordinates. 
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Figure  11.1a.  Rebar  geometry. 


SSS-DFR-89-10119 


11.3  MULTIVARIABLE-FIELD  REPRESENTATION  AND  x*-PERIODICITY 
CONDITION. 


For  the  case  of  bi-directional  steel  reinforcement,  the  x*-periodicity  condition  is 
stated  in  the  form 


u.' 


(2) 


X,- 


2 


(2) 

=  u>  ^ 

I 


fir  A, 


X. 


xi.x* 


2  '  "2'  3J 


.(2) 


JlrA, 


X,  X*, 


'3J 


(2) 


X,  X 


1'  2  '  ^3J 


(11.3) 


u.' 


(2) 


firA, 


X,  Xj,  X*, 


2  J 


(2) 

=  u> 

I 


X»  Xj  ,  Xj; 


fir  A, 
♦  _ 5- 


11.4  TRIAL  DISPLACEMENT  FIELD  (DISPLACEMENT  MICROSTRUCTURE). 

Two  basic  premises  are  now  stated  regarding  the  displacement  microstructure. 
These  are: 


(1)  The  interaction  between  the  two  sets  of  rebar  layouts  which  occupy  the  regions 

a(1) 


and  are  negligible  (this  implies  that  the  interaction  body  force  between  u 
and  is  negligible). 

(2)  The  matrix  displacement  microstructure  can  be  decomposed  as  follows: 


(1) 


uP^(x,x*)  =  U.'^VxCx) 


(2) 

I  (o) 


♦  e 


uj[J](x,xJ,  X*]  .  ,  o[,2]^. 


(1) 


(11.4) 


where  the  superscripts  I  and  II  refer  to  rebar  sets  I  and  II  (see  Fig.  11.1). 

Under  the  assumptions  (1)  and  (2)  above,  the  trial  displacement  fields  for  the 
rebar  and  concrete  are  stated  in  the  form 


Rebar  Set  I  (material  1;  o=l) 

uP(x,x*)  =  Up^x)  ♦  e[?j(x)g^^^^[xj.x*]  ♦  §|(x)h^(^)[x5.xj)]  (11.58) 
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Rebar  Set  II  (material  3;a=3) 

uP^(x,x*]  =  Up^(x)  ♦  e 

Concrete  (material  2;o=2) 

up^[x,x^]  =  Up^(x)  ♦  £ 

[?5(x)^f^’(x5,x*)  .  i?(x)t®(xj,x*] 

where 

i^(x)  =  ij(x) 

(11.6a) 

iP(x)  .  J“(x)  . 

(11.6b) 

A  typical  cubic  cell  is  now  considered  as  illustrated  in  Fig.  11.2.  Further,  a 
multiple  concentric  cylinder’s  approximation  for  this  ceil  is  adapted  as  depicted  in 
Fig.  11.2.  (The  length  of  this  ceil  is  adjusted  to  retain  the  original  volume  of  the 
matrix).  Within  the  context  of  the  concentric  cylinder’s  approximation,  the  trial 
displacements  (11.5)  are  rewritten  in  the  form 

Rebar  Set  I 

up)[x.x*]  =  Up^x)  ♦  e 

[^j(*)g^^^(rj)cos^j  ♦  (rj)sinflj]  ♦  o(e^] 

• 

=  Up^(x)  +  e 

A-«0 
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where 


.(1) 


Rebar  Set  II 


up)[x.x*]  =  up)(x)  ^ 
sindjjj  •♦•  o(e^] 


=  U^^^(x)  ♦  € 


if  (x,^  • 


(3) 


♦  o(£2)  , 


where 


h(3)rr  .  -IlL 
^  ■  ^(3) 


Concrete 


uP>(x,x*)  .  Up>(x)  ♦  «[i?(x)g(^>(rj)cos9j  .  i^i(x)i'>  ('•l)''"‘'l] 
.  e[F(x)h®(rjj)cos9jj  »  if  (x)h'^>(rjj]sin9jj] 

.  o(e2) 


=  (x)  ♦  6 


2l I 

1  1 

f  %  \ 

*  ’‘2 

5i  (x; 

il  -  nWj 

■  *2  *  *2  *2 
*2  *  ’‘a  J 

♦  5T 


r  X* 

*  *3 

ll  -n(«J 

'  *  ^*2  ^  *2 

1  *2  *3  J 

(11.7b) 


(11.8a) 


(11.8b) 
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where 


(^n)  = 


(11.9a) 


(11.9b) 


11.5  PRINCIPLE  OF  VIRTUAL  WORK  FOR  SYNTHESIZED  FIELD. 

In  a  manner  similar  to  that  described  for  a  single  rebar  set,  the  principle  of 
virtual  work  for  the  synthesized  field  can  be  written  as 


a=l  ' 


1 

♦  — 
6 


<IV.  .  7 
(a)  ‘J  'J  *  ' 


a<i) 


-  6u[^^]dSjdx* 


0V. 


(3) 


dV 


I 


a-l 


(a) 


dV 


I 


a=l 


(fl) 


I  I 


(11.10) 


dS 


where  dV  =  dx^dx^dx^,  and  dV,  =  dx*dx*dx*. 


A'tt3 
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and  where 


e.  ■ 
U 


u.  . 

♦  u  •  . 

1 

♦  ~ 

■ 

u 

u 

] 

1  .J 

J , ' 

€ 

L 

J,'  J 

J 

(11.11) 


5u.  .  (5u.  . 

'  iJ  J  / 1 


(11.12) 


11.6  MIXTURE  EQUATIONS  OF  EQUILIBRIUM. 


The  trial  displacement  functions  (11.7)  -  (11.8)  are  now  substituted  into  (11.11) 
and  (11.12)  and  only  those  terms  up  to  an  including  0(1)  are  retained.  The  result  is 
then  substituted  in  (11.10).  together  with  the  displacement  functions  where  required. 
Upon  again  retaining  only  0(1)  terms,  one  thus  obtains  the  following  form  of  the 
principle  of  virtual  work  after  appropriate  integrations  by  parts: 


dS 


f„(l)/la)  ,  „(1)^(13)  ^  pll^(l)  ^  rj3)^(3a)  ^  „(3)^(3a)  ^  pin^(3) 
t  j I  ,j  I  iJ  I  I  j  I  ,j  I  I  J  I 


fn(2)a(2a)  .  „(2),(2a)  .  pi  .  pll)^(2)  _ 

t  J I ,  j  I  I  I  J  I 


,(la)  _  „(2a) 

12  J  _  ^(l)''l2 


*  n(2)Ql  .  n(2)Ql 
n  *<12  *  "  ‘<13 


.(2) 


,(la)  .  _n^ _  ^(2a)  ^  (2)qI  ^  „(2)qI 

22  -  (1)  ‘^11  "  ‘<22  **23 

1  -  n 


dl 


24“)  . 


.  „(2)oI  .  -  2-n.E  .  „(2)|,I  .  „(2)„I 

"  *<33  1  (1)  ^^23  "  *^23  "  ‘<32 

1  -  n  ' 


(4  =  4) 
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(la)  _  _ .(2a)  ^  (2)^1  ^  (2)„I 


'31 


t  (1)  31 

1  -  n '  ' 

.(2) 


13 


'12 


_ _(2a)  (2)f,I  ^  -(2)rI 


(la)  _nj 

'33  -  j  _  ^(1)“33 


33 


32 


''^3 


.P) 


2  .  „P)pn  .  ^(23)  .„(2,pii  ,  „(3,,n 


=  aiH] 


'  (33)  _  ji?L_  423)  .  42)41 . 42)41 


'23 


1  -  3(5)  23 


'21 


(33)  _  .jP)  (2a)  (2)  .11  .  (2).II  JII 

(3)  "33  '33  '31  °^3 


'33 


1  -  n' 


ill 


,(3a) 

^11 


,(2) 


1  -  n(3)  ‘'ll 


^(2a)  (2) oil  (2) oil 


.(2) 


(3a)  ^  (2)rII  ^  _n  _ 

^12  ♦  n  K23  ^  ^  ^^3j  vj2 


(2a)  ^  (2)pII 
m  ▼  n  1^2  j 


3jn 


dV  s  0  .  (11.13) 


The  Euler-Lagrange  equations  of  (11.13)  lead  immediately  to  the  following 
mixture  equilibrium  equations  and  associated  internal  condensation  relations  which  can 
be  used  to  eliminate  the  variables  and  S^. 


Mixture  Equilibrium  Equations. 

nd),!?*?  .  .  P?  =  0  , 

J‘,J  »  < 

n(3)a(3»)  .  n^^if^^a)  ,  pH  =  q  . 
J'.J  '  « 

n(2)^(2a)  ,  „(2),(2a)  .  pi  .  pll  ,  0 
J'.J  '  • 


(11.14a) 

(11.14b) 

(11.14c) 


Internal  Condensation  Relations. 
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^(la) 

(2a)  , 

"‘"’K  ♦  <>13  ] 

—  n 

(11.15a) 

^12  ' 

l-n(l) 

-  y  I 

^(la) 

a(2a^  . 

a  (3)  (4  •  4) 

=  0  . 

(11.15b) 

^22  ■ 

1  -  n(l) 

(2a)  , 

[rI  ^  qI 

*  ^23  *  *^32]  "  °  ' 

(11.15c) 

*^23  ■ 

1  -  23 

2  r22  ^33 

n(2)  (2a) 

.  ti®[Rl3  .  rJj] 

1=0. 

(11.15d) 

""31 

i.„a)‘'3i 

^(la) 

t.<^>  (2a) 

A 

(11.15e) 

^33  ■ 

l.n(l)  33 

-  U  > 

.(3a) 

0^2)  (2a) 

E®  foil  .  rII 

nil  nlll  A 

■"  *^31  *^13)  *  ®  ' 

ril  IRal 

^^31  - 

l.n(3)"31 

2  ril  *^33 

X . XDa  J 

^(3a) 

n^2)  (2a) 

—  n 

(11.16b) 

^^23  ■ 

l.n(3)"23 

•  y  / 

a  (3**) 

n^2)  (2a) 

-  *  oJl] 

—  n 

(11.16c) 

^^33 

■  1  -  oW  33 

1  -  u  , 

.(3a) 

n®  (2a)  , 

(11.16d) 

"'ll  • 

1-n®  11 

-  u  , 

^(3a) 

ti®  (2a) 

«  n 

(11.16e) 

^12  - 

l.n(3)"l2 

-  u  , 

Boundary  Conditions. 

«U<*)  =  0  or  (11.17) 

I  J  I  • 


A-Bb 
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Averages 


(fM  =  - 1 - 

"ij  -„(«),3/2 


,(a) 


ij  * 


(11.18a) 


[rI  Rn]  _  ■_  1  __ 

rij'^ijJ  -  „(2)^3/2 


(2) 


L  i,j*'  ij  * 


(11.18b) 


°  „(2),3/2  W,j"  ' 

"  '  JJJu(2) 


(11.18c) 


pl.  » 


i  -  „3/2 


t?*dS* 


(1) 


pii _ 


dS* 

I 


av 


(3) 


(11.18d) 


(ll.lSe) 


in  which  dS*  is  an  infinitesimal  interface  area  element  on  or  avj^^’  and 


ff®^dV,  ,  (11.18f) 

v(«)' 

Yp^dV^  (ll.lSg) 

M 

11.7  MIXTURE  CONSTITUTIVE  RELATIONS. 

•  In  order  to  generate  appropriate  mixture  constitutive  relations  to  accompany  the 

equilibrium  equations  (11.14),  one  can  utilize  a  mixed  weighted  residual  procedure.  The 
situation  here  is  similar  to  that  encountered  for  the  case  of  a  single  rebar  set. 


Y(m)  _  1 

'  -  „C),3/2 


f(M) _ 1 

T  •  — 


A-87 


SSS-DFR-89-10119 


Trial  Displacement  Field. 

The  trial  displacement  field  (11.7)  -  (11.9)  is  again  used  in  the  mixed  residual 
procedure.  This  field  can  be  expressed  in  the  form 


u^“^(x,x»;e)  =  U^“^(x)  ♦  eG^“^(x*)S 


(11.19a) 


where 


This  trial  field  leads  to  the  following  strain  field  via  (11.11)  when  terms  of  0(1)  only 
are  retained: 


s  e(u^“^]  ♦  Q^®^(x*)S(x) 


(11.20a) 


where 


EtuW)  s  K,i,Uj_j,U3  3,Uj  3  .  U3_3,U,  3  .  U3_,,Uj_j  .  Uj  J (11 . 20b) 
Trial  Stress  Field. 

The  trial  stress  field  to  be  used  in  the  mixed  residual  procedure  is  selected  in 
the  form: 


=  rf®)(x)  >  5„2^i jkUx*)t (*)  *  eaff, 


U 


u 


'kZ 


ij(l) 


(11.21a) 


where 


(( 


)ji 

V 

and  where 


^ijkt(x*)dV.  =  0  , 

v(2) 

# 


(11.21b) 


ii(l) 


are  defined  by 
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‘'ll(l) 

(a) 

\  §{x)g^^(“)(x*)  .  1  Pj^x)  .  h^^(“)(x*) 

CM 

CM 

I  -  \  P3(x»)h^(“)(x*) 

^33(1) 

= 

\  P2(*)g^^“^(x*)  -  I  P3(x)h^^“^(x*)  ♦  I  P3^(x)g^^(“^x*)  ♦  JpJ^x) fS^) 

""23(1) 

\  P^(x)h^(“)(x*)  *  \  P3(x)g^(“)(x*)  *  \  P”(x)g^^(“^(x*) 

‘^31(1) 

i  P{(x)^(“)(x*)  .  \  P^^x)^^(“)(x*)  .  \  Pj^(x)J^(“)(x*) 

‘'12(1). 

^J(x)g^“hx*)  >  \  P^^(x)S^(®>  (x*) 

or  by 


'S  °  •  HjW  (.•)?”(») 


where 


p'  ^  [pI.S'.pJ]"  .  p"  -  lp{'.p^'.pjf  ■ 

Thus,  (11.21a)  can  be  rewritten  as 

ff^“^(x,x*)  =  t(“^(x)  ♦  5g2^^^^(x*)t^^^(x)  ♦  eH^J^(x*)P^(x) 


(11.21c) 

(11.21d) 

(11.21e) 


♦eH^5^x*)P^^(x)  . 


(11.22) 


Several  remarks  are  in  order  at  this  point.  First,  the  first  two  terms  of  the  trial 
stress  field  (11.22),  i.e. 


=  T^“^(x)  ♦  T^^^(x*)b^^^6 


was  constructed  to  satisfy  j*~  ®  obtained  from  the  displacement 

in  (11.19a).  Second,  it  is  noted  that 


a2 
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T^^^(x*)dV  =  0  , 


r(2) 


=  0 


,(a) 


=  0 


SO  that 


,(2) 


i.e.,  is  the  average  stress 

Finally,  the  forms  and  introduced  where  constructed  to  satisfy  the 
integrability  condition  for 


3/2 


8V, 


(1) 


^(1)  ha  -  _ 


3V. 


(3) 


.#dA„  s 

jt(i),j*  II  t 


Local  Constitutive  Relations. 

The  local  constitutive  relations  for  each  component  material  are  specified  in  the 
rate-form 


(a);  (a) 


(11.23) 


where  is  a  6x1  matrix,  is  a  6x6  matrix,  and  is  a  6xi  matrix.  The 
coefficients  are  tangaui  moduli.  Thus.  (11.23)  incorporates  incremental  plasticity. 
All  stresses  in  (11.23)  are  transverse  stresses. 

Mixed  Weighted  Residual  Form. 


The  appropriate  form  of  the  mixed  weighted  residual  for  the  present  problem  is 


C 


dV 


4 
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rr 


1 

♦  — 
e 


(2) 


""ds’dx*  ^  ^ 


II  j  *  1 


6v, 


(1) 


ai 


(2)  . 


5*(3)jf*II^S*jdx2 


6V, 


(3) 


I 


a=l  j 


5;(“)Tf(“)dV 


rW 


dV  ♦ 


r  r  2 


6v, 


I 


l  0=1  J 


(a) 


dS, 


where 


L(a)  , 


® 11 ' ®22 ' *33 ' ^®23 ' ^®31 ' ^® 12 J 


1(0)^ 


(11.24) 


(11-25) 


Homogenized  Rate  Constitutive  Relations. 

Upon  substitution  of  the  fields 

;(«)  ,0(«)(.,  .  ee<">(.')5(x)  . 

;(“)  =  e(0(“>)  .  R<“)  («*)5  ,  (11-26) 

J(«)  ^  f(«)  »  .  eH{“\x*)^(x)  .  cH|{>(x*)f^^(x) 

into  (11.24),  one  obtains,  as  the  Euler-Lagrange  equations,  the  rate  form  of  (11.14)  - 
(11.16)  together  with  the  rate  constitutive  relations 


v; 


Q(®)(x*)dV 

(a) 


r  r  r  r 


§- 


i_ 

V* 


C^“^dV 


,(<*) 


;.(a) 


-  6 


02 


;(a)T(2)dV„ 


(a) 


t(2)  .  e 


c(«)H(®)dV^ 
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ABSTRACT 


A  new  internal  variable  theory  of  plastic-fracturing  solids  is  described.  The  key 
concept  upon  which  the  theory  rests  is  a  mapping  that  transforms  the  current, 
damaged  and  generally  anisotropic  state  of  a  material  into  an  isotropic  undamaged 
state.  The  effect  of  cracking  on  the  constitutive  behavior  is  accounted  for  through  a 
second-order  integrity  tensor  which  is  governed  by  an  evolution  equation.  The 
integrity  tensor  also  allows  anisotropy  to  be  described  in  a  convenient  way  without 
introducing  a  fourth-order  tensor,  as  is  the  usual  approach.  The  model  describes  both 
stiffness  degradation  and  yield  limit  degradation  due  to  developing  damage,  as  well  as 
dilatancy.  and  contains,  as  special  cases,  several  models  that  are  noteworthy,  including 
the  classical  elastic-fracturing  model  and  the  endochronic  plasticity  model.  Finally,  the 
model  is  applied  to  recent  data  on  the  response  of  plain  concrete  to  simple  tension  in 
the  presence  of  developing  damage. 
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INTRODUCTION 


Many  of  the  important  questions  faced  by  the  Air  Force  today  regarding  the 
hardness  of  concrete  defense  structures  require  a  knowledge  and  understanding  of  the 
likelihood  that  a  given  structure  will  fail  under  a  prescribed  loading  environment.  To 
address  this  issue  in  the  most  cost-effective  manner,  analytic  methods  must  be 
available  to  allow  system  designers  to  systematically  vary  design  parameters  in  an 
expedient  and  efficient  manner  so  that  an  optimum  structural  design  can  be  achieved. 
Finite  element  methods  are  now  available  for  this  purpose,  but  the  results  obtained 
with  these  methods  can  only  be  as  accurate  as  the  constitutive  relations  used  in 
conjunction  with  them.  To  examine  the  hardness  of  concrete  structure  requires  a 
constitutive  model  of  plain  concrete  which  not  only  is  capable  of  describing  the 
inelastic  behavior  of  intact,  unfractured  material  but  is  capable  of  describing  developing 
damage  and  cracking  and  their  effect  on  the  overall  constitutive  behavior.  For  a 
constitutive  model  to  be  useful  in  addressing  these  important  issues,  it  must  be  able 
to  treat  both  elastic,  plastic  and  cracking  behavior. 

Recently,  Valanis  and  Read  (1986)  developed  an  advanced  nonlinear  constitutive 
model  for  plain  concrete,  on  the  basis  of  the  endochronic  theory  of  plasticity.  The 
model  was  applied  to  an  extensive  set  of  laboratory  data  for  a  medium  strength  plain 
concrete  reported  by  Scavuzzo,  et  aL  (1983).  The  test  program  consisted  of  a  variety 
of  complex  loading  paths  designed  to  explore  different  facets  of  concrete  behavior. 
The  stress  paths  were  selected,  however,  so  that  no  significant  (macroscopic)  cracking 
of  the  concrete  occurred  during  the  tests.  As  shown  by  Valanis  and  Read  (1986), 
the  endochronic  model  for  concrete  was  remarkably  successful  in  describing  the 
behavior  of  the  concrete  over  the  wide  variety  of  complex  stress  paths  explored  in  the 
laboratory. 

The  constitutive  model  for  concrete  mentioned  above  has  several  limitations, 
however,  which  restrict  its  application  to  stress  states  below  failure.  First,  the  model 
is  incapable  of  accounting  for  the  effect  of  macrocracking  on  material  behavior. 
Secondly,  since  the  model  is  isotropic,  it  is  unable  to  account  for  the  anisotropy  that 
develops  when  significant  macrocracking  occurs  in  preferred  directions.  Finally,  in  the 
case  of  shearing  at  fixed  hydrostatic  pressure,  the  model  exhibits  only  compaction;  the 
dilatant  behavior  observed  when  cracked  concrete  is  sheared  cannot  be  described  by 
this  model. 


Recently.  Valanis  (1987,  1988a)  developed  an  approach  for  extending  the 
endochronic  theory  to  include  the  effects  of  damage  and  fracture.  The  result  is  a 
new  endochronic  plastic-damaging  model  which  possesses  all  of  the  attractive  features 
of  the  earlier  concrete  model  (Valanis  and  Read,  1986)  but  which  also  can  treat 
cracking,  anisotropy  and  dilatancy.  The  resulting  model  appears  to  have  the  desired 
features  for  describing  the  complete  spectrum  of  concrete  response.  Including  plastic 
flow,  cracking  and  dilatancy.  Anisotropy  is  accounted  for  in  a  clever  way  through  the 
use  of  the  second-order  integrity  tensor  rather  than  by  a  fourth-order  tensor,  as  is 
the  usual  approach.  The  use  of  a  second-order  tensor  for  this  purpose  greatly 
reduces  the  analytic  complexity  of  the  model  and  simplifies  the  experimental 
determination  of  the  material  parameters.  The  model  describes  both  stiffness 
degradation  and  yield  limit  degradation  due  to  developing  damage  and  contains,  as 
special  cases,  several  models  that  are  noteworthy,  such  as  the  classical  elastic¬ 
fracturing  model. 

Parenthetically,  it  should  be  noted  that  since  the  endochronic  theory  contains 
classical  plasticity  as  a  special  case  (as  shown  by  Valanis.  1980),  the  approach 
described  herein  can  also  be  used  to  determine  the  form  taken  by  classical  plasticity 
theory  when  there  is  cracking,  although  this  has  not  been  done  in  the  work  reported 
herein. 

The  present  study  represents  the  first  attempt  to  explore  the  application  of  the 
new  endochronic  plastic-fracturing  theory  to  plain  concrete  and  to  develop  a  numerical 
approach  for  dealing  with  the  system  of  governing  equations  computationally.  After 
the  basic  formulation  of  the  new  theory  is  given,  attention  is  focussed  on  the  special 
case  of  uniaxial  loading  to  elucidate  some  of  the  basic  features  of  the  theory.  The 
model  is  applied  to  the  problem  of  the  response  of  plain  concrete  to  uniaxial  tension 
in  the  presence  of  developing  damage,  using  the  recent  laboratory  data  reported  by 
Gopalaratnam  and  Shah  (1985).  It  is  shown  that  the  resulting  model  describes  the 
observed  behavior  very  well,  including  the  post-peak  response  and  the  unloading¬ 
reloading  behavior.  The  observed  hysteresis  during  the  unload-reload  cycles  is, 
however,  somewhat  greater  than  the  model  predicts. 
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Section  2 


FORMULATION  OF  MODEL 


In  this  section,  the  basic  formulation  of  the  endochronic  plastic  fracturing  model 
is  given,  following  the  work  of  Valanis  (1988a).  The  key  concept,  upon  which  the 
formulation  rests,  is  a  mapping  that  transforms  the  current,  damaged  and  generally 
anisotropic  state  of  a  material  into  an  isotropic  undamaged  state,  which  is  called  the 
"transformed  state".  Since  the  standard  endochronic  plasticity  theory  applies  to 
isotropic,  undamaged  material,  it  is^taken  here  as  the  constitutive  model  for  the 
material  in  the  transformed  state.  Then,  by  using  the  mapping  from  the 

transformed  state  to  the  current  state,  the  constitutive  relations  for  an  endochronic 
plastic-fracturing  model  are  found.  The  effect  of  cracking  is  reflected  through  a 
second-order  integrity  tensor  jj,  which  is  governed  by  an  evolution  equation. 

2.1  MAPPING  FROM  CURRENT  STATE  TO  TRANSFORMED  STATE. 


To  begin,  we  note  that  damage  and  fracture  produce  a  reduction  in  material 
integrity.  With  this  in  mind,  an  integrity  tensor.  jJ,  is  introduced  which  has  the 
following  properties:  it  is  a  second-order  symmetric  tensor  which  is  equal  to  the  unit 
tensor  jj  in  the  undamaged  virgin  material  and  to  the  null  tensor  in  fully  fractured 
materials,  i.e..  material  which  cannot  support  tensile  stress  in  any  direction.  The 
integrity  tensor  is  positive  semi-definite  but  can  be  regarded  as  positive  definite  in 
the  sense  that  its  norm  may  be  made  as  close  to  zero  as  one  wishes,  Vi^ithout  it 
actually  being  zero;  this  proves  to  be  a  useful  concept  for  the  purposes  of  the  ensuing 
analysis.  Inasmuch  as  is  symmetric  and  may  be  regarded  as  positive  definite,  it  can 
be  represented  in  the  polar  form 


<  =  B 


where 


exists  and  is  unique,  i.e.. 


(2-1) 


In  indicial  notation,  Eq.  (2-1)  reads 


'J 


(2-2) 


(2-3) 


***  Other  constitutive  models,  such  as  those  from  classical  plasticity,  for  example, 
could  also  be  used  to  describe  the  material  in  the  transformed  state. 
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On  the  basis  of  the  above  considerations,  the  free  energy  density  ^  for  a  plastic¬ 
fracturing  solid  may  be  expressed  in  the  following  form: 


~  2  ^  ^iike.(^)Kg.  ■  ■  ‘^ij) 

r 

where 

=  S'  ^ij  (2-5) 


Here,  in  the  formalism  of  internal  variable  theory,  the  integrity  tensor  is  an  internal 
variable  which  characterizes  damage,  while  the  internal  variables  characterize  plastic 
deformation. 

The  stress  g  is  obtained  from  ^  in  the  usual  manner,  i.e.. 


which,  on  the  basis  of  £q.  (2-4).  leads  to  the  result: 

‘^ij  “  X  {^iSj^ke.  ^2  ^ikS£}Ke  "  \e} 


(2-7) 


We  now  consider  a  mapping  £  =  £(x)  from  the  current  damaged  state  to  the 
"transformed"  state,  as  shown  in  Figure  1.  If  g  is  taken  as  the  transpose  of  the 
deformation  gradient  £  between  the  current  and  transformed  states,  i.e.. 


B  =  (2-8) 

then  the  strain  tensor  and  internal  variables  in  the  transformed  state  can  be 
defined  in  terms  of  the  strain  tensor  g  and  internal  variables  in  the  current  state 
as  follows: 


e.o  =  B.  . 

k£.  ki  £.j  ij 


(2-9) 


or 
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Transformed  Undamaged  State 


(2-10) 


and 


q^o  =  B.  -Bo .  qT . 

^k£.  k  I  K.j  ^  I  j 


or 


‘’ki  ■  ®ki®lj  ‘’ij 


(2-11) 


(2-12) 


The  use  of  Eqs.  (2-2),  (2-3)  and  (2-9)  to  (2-12)  in  Eq.  (2-4)  "normalizes"  this  last 
equation  and  leads  to  the  following  simple  expression: 


kZ 


'’kJ 


(2-13) 


which  implies  that  the  material  in  the  transformed  state  is  isotropic  and  undamaged. 

The  stress  g  in  the  transformed  state  may  be  obtained  by  differentiating 
Eq.  (2-13)  with  respect  to  leading  to  the  expression: 


"ij  '  X]{Vij[^k  -  Okk)  *  *2hj  -  ‘'ij)} 


Upon  comparing  Eqs.  (2-7)  and  (2-14).  we  find  that 


a. 0  =  BT?  BTJ  a. . 
kZ  ik  jt  IJ 


or 


'ke. 


=  B. ,  B.<,a. . 

ik  j£.  Ij 


(2-14) 

(2-15) 

(2-16) 
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In  view  of  the  above  developments,  it  is  concluded  that  the  thermodynamic  analysis 
may  be  done  in  the  transformed  state,  where  the  deformed  material  is  both  isotropic 
and  undamaged,  and  then  mapped  back  to  the  current  deformed  state,  in  which  the 
material  is  both  damaged  and  anisotropic,  by  means  of  the  transformations  (2-9). 
•  (2-10).  (2-15)  and  (2-16). 

2.2  CONSTITUTIVE  EQUATION  FOR  THE  TRANSFORMED  STATE. 

Consider  now  the  evolution  equation  for  the  plastic  internal  variables  in  the 
^  transformed  state,  which  we  take  In  the  usual  form: 


b^ 

ag’’  ~ 


0 


(2-17) 


Here,  b*^  is  an  isotropic  tensor  of  the  fourth  order  (the  ‘’resistance’’  tensor),  in  which 
case  it  has  the  following  representation: 


(2-18) 


r  ^ '  r 

Note  that  the  tensor  b  in  the  current  state  is  related  to  b  in  the  transformed  state 
by  the  tensor  transformation: 


bC..o  =  B.  B.  B.  .Bo  b'"  . 
ijke.  ir  js  kt  £x:^  rstq 

Therefore,  in  view  of  Eqs.  (2-18)  and  (2-19),  it  follows  that: 


(2-19) 


^ijk£.  ”  ^l^ij^kt  ^  ^2^ik^j£. 


(2-20) 


which  shows  that  the  integrity  tensor  ^  affects  the  resistance  tensor  b^ 

The  solution  of  Eq,  (2-17)  for  a  generic  r  is  given  elsewhere  (see  Valanis  and 
Read,  1989),  and  will  not  be  repeated  here.  Briefly,  is  found  from  Eq.  (2-17), 

using  Eqs.  (2-13)  and  (2-18).  The  result  is  substituted  into  Eq.  (2-14),  leading  to 
the  following  result: 
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(2-21) 


where 


(2-22) 


=  ll‘'s’’||^  ♦  ‘‘^hkk)' 


(2-23) 


Here,  F-  and  F^  are  hardening  functions  for  deviatoric  and  hydrostatic  response  which 
satisfy  tne  conditions: 


Fn  >0  I  Fy  >  0 


(2-24) 


Furthermore,  the  double  bars  around  a  symbol  denote  its  Euclidean  norm,  and 


-“r"D 


(2-25) 


KfzJ  =  >  K_e 


where  /i^,  K^,  and  are  ail  positive  constants.  Specifically: 


^  =  *2 


=  Aj  .  i 


^  •  5  ‘■2] 


(2-26) 
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The  above  constitutive  equation  can  alternately  be  expressed  in  terms  of  plastic 
strain  integrals  with  singular  kernels.  The  details  for  accomplishing  this  are  also  given 
by  Valanis  and  Read  (1989).  Thus,  Eq.  (2-20)  may  be  written  in  the  form: 


o 


J 

o 


(2-27) 


where  p(z)  and  ^(z)  are  weakly  singular  kernel  functions  that  satisfy  the  conditions 
p{0)  =  ^(0)  =  «  such  that  the  integrals 


’^D 

p  (y)  dy 


^(y)dy 


exist  for  all  Zq  ^  0  and  \  0. 

To  complete  the  specification  of  the  constitutive  model  for  the  transformed  state, 
we  must  add  Hooke’s  law: 


■kkJ^ij 


(2-28) 


where  X  and  p  are  the  Lame  constants. 

2.3  CONSTITUTIVE  EQUATION  FOR  THE  CURRENT  STATE. 

The  relation  between  g  and  g  in  the  current  deformed,  damaged  configuration  is 
found  by  using  Eqs.  (2-8)  and  (2-16)  in  Eq.  (2-26),  with  the  result: 


a  =  B  .B  . 
rs  n  sj 


'’(^0  -  'IdT 


PH 


rs 


(2-29) 


O 


8-14 


To  eliminate  the  tensor  g  from  the  above  expression  in  favor  of  the  tensor  ^  we 
proceed  according  to  Valanis  (1988b)  in  the  following  manner.  Since  is  a 
symmetric,  second-order  tensor,  it  can  be  expressed  in  the  form: 


=  N 

Jw  ^  JCr  M 


(2-30) 


where  N  denotes  the  matrix  of  the  eigenvectors  of  and  is  the  (diagonal)  matrix 
of  the  eigenvalues  of  jj.  Then,  in  view  of  Eqs  (2-1)  and  (2-30).  g  can  be  expressed 
in  terms  of  ^  as  follows: 

B  =  (2-31) 

which,  in  indicial  notation,  reads 

Bij  =  Kk)"'"  N-k  P-32) 


Thus,  on  the  basis  of  Eq  (2-32),  Eq  (2-29)  can  be  expressed  solely  in  terms  of 

It  should  also  be  noted  that  there  are  conditions  under  which  Eq  (2-29)  reduces 
to  a  more  simplified  form.  Consider,  for  instance,  the  following  conditions: 

(i)  §  is  a  slowly  varying  function  of  z 

(ii)  p  and  ^  are  rapidly  decaying  functions  such  that  F(z')  may  be 
evaluated  at  z'  =  z  in  the  integrals  on  the  righthand  side  of 
Eq.  (2-29). 

In  view  of  (i),  we  can  write: 


dz'[®ik®j£.®k< 

B  B 

J  ix  jt  dz' 

dz'  lhe.^kd 

-  dz' 
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Therefore,  on  the  basis  of  Eqs.  (2-33).  Eq.  (2-29)  can  be  written  as 


=  - 


(2-34) 


KsKil  ^("h  - 


dz 


Consider  now  the  expression  (2-23)  for  the  intrinsic  time  scale  d^.  Upon  using 
Eqs.  (2-8)  and  (2-9).  Eq.  (2-23)  transforms  to  the  following  expression  for  df  in  the 
current  state: 


df^s  ^  ^  de^  deP  ♦  k^  de^  1‘ 

^  ^mp  ^qn  pq  mn  r  mn  mnJ 


(2-35) 


Note  that  the  tensor  de**  in  the  current  configuration  is  not  deviatoric.  i.e., 
tr(dg*’)  ^  0:  this  follows  from  the  fact  that  the  current  state  is,  in  general, 
anisotropic.  To  see  this,  recall  the  expression  for  the  deviatoric  plastic  strain  g*’  in 
the  transformed  state: 


's'"-?  -I 


(2-36) 


Then  from  Eqs.  (2-9)  and  (2-10).  we  can  write,  using  matrix  notation: 


(2-37) 


Substituting  Eqs.  (2-37)  into  (2-36)  leads  to  the  result: 


'  s”  -  5  B'Mb'T 


(2-38) 
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It  can  be  shown  that 


tr[l^]  =  tr[i  eP]  (2-39) 

Thus.  Eq.  (2-38)  can  be  written  as 

-  I  tr(i  gP]  (2-40) 

since,  from  Eq.  (2-1).  we  have  =  §  ^(§  As  a  result,  it  follows  that 

tr[eP]  =  tr[£P]  -  |  tr[jj  £p)tr(^*^]  (2-41) 

An  inspection  of  Eq.  (2-41)  shows  that  tr(gP)  =  0  only  when  [J  =  g,  i.e.,  in  the 
initial  undamaged  (isotropic)  state. 

To  express  the  righthand  side  of  Eq.  (2-35)  solely  in  terms  of  the  tensor  dgP, 
we  proceed  as  follows.  First,  it  can  be  shown  from  Eq.  (2-40)  that 

<  •  s”  =  <  •  -is"  -  5  <>'•(<  •  ■i£’’)iS  (2-«) 

Equation  (2-35)  can  be  written  in  the  following  matrix  form: 

-  tr{(<  .  ds^JIds”  .  .  df"]}^  (2-43) 

Upon  substituting  Eq.  (2-42)  into  (2-43),  we  find: 

ds^  =  tr{(<  .  dj'’)(d£'’  •  .  (k^  -  l]{tr(<  .  djPjj^  (2-44) 

In  the  isotropic,  undamaged  transformed  configuration,  Hooke's  law  has  the 
following  form,  which  was  given  earlier  in  Eq.  (2-28): 

2  =  X  tr(g  -  gP)g  ♦  2/i[g  -  gP)  (2-45) 


This  expression  was  first  developed  by  H.  Murakami. 
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To  express  Eq.  (2-45)  in  the  current  damaged  state,  we  introduce  the  transformations 
given  by  Eqs.  (2-9)  and  (2-16).  Proceeding  in  this  manner,  it  follows  that 

£  =  X  i  tr[t[£  -  jP]]  .  2,1  i(£  -  £'’]i  (2-46) 

which,  in  indicial  notation,  reads 

^ij  "  ^  ^k^jiKe  ■  ^kJ 

For  many  geomaterials,  including  concrete,  the  hardening  functions  are  typically 
of  the  general  form: 

fH  = 

(2-48) 

=  foH 

where  £**  and  a  denote,  respectively,  the  plastic  volume  strain  and  the  hydrostatic 
pressure  in  the  current  configuration. 

To  summarize  the  main  results  from  this  section,  the  system  of  equations  which 
govern  the  behavior  of  the  endochronic  plastic-fracturing  model  in  the  current  damaged 
state  consist  of  Eqs.  (2-22),  (2-29),  (2-44),  (2-46)  and  an  evolution  equation  for 

The  key  to  the  success  of  the  damage  model  described  above  lies  in  having  the 
appropriate  evolution  equation  for  ^  for  the  material  of  interest.  Ideally,  it  would  be 
desirable  to  have  the  evolution  equation  reflect  the  underlying  micromechanical  fracture 
mechanisms.  A  research  effort  to  develop  such  an  evolution  equation  for  ^  is 
currently  being  conducted  by  Valanis  in  a  separate  AFOSR-funded  program,  but  the 
results  from  this  study  are  not  expected  to  be  available  for  sometime. 

In  the  following  section,  the  response  of  the  model  to  the  special  case  of 
uniaxial  loading  is  considered  in  detail  to  explore  some  of  the  general  features  of  the 
theory.  The  model  is  then  applied  to  the  recent  data  by  Gopalaratnam  and  Shah 
(1985)  on  the  response  of  plain  concrete  to  simple  tension  in  the  presence  of 
developing  damage  and  unloading-reloading  cycles.  For  this  purpose,  an  evolution 
equation  suggested  by  the  data  is  adopted  and  successfully  used. 
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Section  3 


UNIAXIAL  LOADING 

In  this  section,  the  case  of  uniaxial  loading  is  considered  in  detail  to  illustrate 
some  basic  features  of  the  endochronic  plastic-fracturing  model  described  in  the 
preceding  section.  First,  the  governing  equations  for  uniaxial  loading  are  deduced  from 
the  general  equations  given  in  Section  2.  Then,  several  interesting  features  of  the 
model  are  discussed,  including  the  several  important  special  cases  that  the  model 
reduces  to  under  certain  conditions.  Following  this,  an  incremental  scheme  is 
developed  for  integrating  the  system  of  governing  equations.  Finally,  the  model  is 
applied  to  the  recent  laboratory  data  of  Gopalaratnam  and  Shah  (1985)  for  the 
response  of  plain  concrete  to  uniaxial  tension,  and  it  is  shown  that  the  model 
provides  an  excellent  description  of  the  observed  behavior. 

3.1  GOVERNING  EQUATIONS. 

Consider  the  case  of  uniaxial  tension  in  which  the  applied  stress  is  directed 
along  the  x.-axis.  In  this  case,  we  have  B..  =  0  for  i  /  j  and  a  .  =  a  From 
Eqs.  (2-9)  and  (2-16)  we  therefore  find  that  * 

A 

‘^ll  *  ^11  ^11 

hi  =  ^11  hi  (3-1) 

eP  =  ^  eP 

^11  Pll  ^11 


The  constitutive  behavior  of  the  endochronic  plastic-fracturing  model  in  the 
transformed  state  is  described  by  the  system  of  equations  given  in  Section  2.2  which, 
for  the  case  of  uniaxial  loading  and  plastic  incompressibility,  reduce  to  the  following 
form: 


hi  = 


E(2 


deP 


(3-2) 
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^11  "  ■  ^ll] 


(3-3) 


= 


3 

2 


(3-4) 


dz  =  ^  (3-5) 


Here.  E(z)  is  a  weakly  singular  kernel  function.  denotes  Young's  modulus,  z 

is  the  intrinsic  time  and  ^  is  the  intrinsic  time  scale.  In  addition.  F  represents  the 
shear  softening  function,  and  the  vertical  bars  surrounding  a  symbol  denote  its 
absolute  value.  Note  that  Eqs.  (3-2)  to  (3-5)  are  the  governing  equations  for  the 
standard  (undamaged]  endochronic  theory. 

The  kernel  function  E(z]  is  given,  in  the  usual  manner,  by  the  Dirichlet  series: 


(3-6) 


where,  in  order  to  satisfy  the  conditions  imposed  by  the  weakly  singular  nature  of 
E(z),  i.e..  E(0)  =  •  and  E(y)dy  <  0,  the  positive  constants  and  must 
satisfy  the  following  conditions: 

m 

r=l 

(3-7) 


***** 


As  shown  by  Valanis  and  Read  (1989).  E(z)  is  related  to  p(z)  according  to 
the  expression  E(z)  =  3/2  ^(z). 
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From  experience,  we  have  found  that,  for  most  materials,  only  two  or  three  terms  of 
the  series  (3-6)  are  necessary  to  provide  an  adequate  modeling  capability.  In  view  of 
this,  we  write 


E(2)  =  ^ 

r=l 


z 


(3-8) 


where  n  is  finite  and  care  is  taken  to  insure  that  the  infinitely  large  value  of  E(0)  is 
suitably  approximated  by  a  large  finite  value. 


The  above  system  of  equations,  which  apply  to  the  transformed  state,  are  now 
mapped  to  the  current  (damaged)  state  through  the  use  of  Eqs.  (3-1).  The  following 
set  of  equations  are  then  obtained  for  the  current  state: 


/Z 

^  E(z  -  z  )  —  dz 

o 

(3-9) 

O’  =  (c  -  e^) 

'  0 

(3-10) 

d^  = , 

(3-11) 

(3-12) 

where  the  subscripts  have  been  suppressed  for  convenience.  To  the  above  equations, 
we  must  add  an  expression  which  defines  the  manner  in  which  ^  depends  upon  the 
deformation. 

In  uniaxial  loading,  it  is  reasonable  to  assume  that  f  will  depend  upon  a  variable 
related  to  the  axial  extension,  and  the  three  likely  candidates  are  the  total  axial  strain 
6.  the  elastic  axial  strain  €*  and  the  plastic  axial  strain  e^.  During  the  course  of 
applying  the  model  to  the  uniaxial  concrete  data  of  Gopalaratnam  and  Shah  (1985) 
(see  Section  3.4),  we  investigated  the  dependence  of  ^  on  each  of  the  three  strain 
variables  noted  above  and  found  that  a  dependence  on  the  axial  plastic  strain  was 
preferred.  The  reasons  for  this  are  as  follows.  For  the  case  in  which  ^  is  taken  to 
depend  on  the  total  strain  £,  it  was  found  that  the  model  was  overly  sensitive  to  the 
form  adopted  for  ^(c).  Very  slight  deviations  from  a  form  which  gave  good 
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agreement  with  the  Gopalaratnam-Shah  data  lead  to  difficulties  in  numerically 
integrating  the  system  of  governing  equations;  this  was  not  the  case  when  ^  was 
taken  to  depend  on  the  axial  plastic  strain,  Furthermore,  on  the  basis  of  the 
results  described  in  Section  3.4,  it  was  found  that  ^  could  not  be  described  uniquely 
as  a  function  of  the  axial  elastic  strain,  e*.  i.e.,  for  a  range  of  values  of  6^,  ^  was 
double-valued.  Again,  this  was  not  the  case  when  ^  depended  on  the  axial  plastic 
strain,  For  the  reasons  cited  above.  ^  is  taken  to  depend  on  e**.  We  therefore 
adopt  the  general  form: 

f  =  C('S)  (3-13) 

where  e**  denotes  the  plastic  component  of  the  total  tensile  strain  e,  i.e.,  e**  =  e  - 
£*.  and  the  subscript  m  indicates  the  maximum  value  attained. 

The  corresponding  evolution  equation  for  ^  is  then 
H(e^deP  ,  if  dfP  >  0  on  ^ 

d^  =  (3-14) 

0  ,  otherwise 

where  H  =  d^/de**.  Thus,  in  the  event  of  unloading  from  a  tensile  strain  state,  we 
have  de^  =  0  and  therefore  d4  =  0,  as  should  be. 

3.2  COMMENTS  ON  THE  MODEL. 

The  model  described  above  contains,  as  special  cases,  several  models  which  are 
particularly  noteworthy,  whose  responses  are  depicted  in  Figure  2.  We  note,  first  of 
all.  that  there  are  ideally  two  kinds  of  damage  (or  degradation)  that  are  observed  in 
materials  (see  Bazant  and  Lin.  1988): 

1.  Degradation  of  material  stiffness,  which  is  caused  by  microcracking, 
nucleation  and  void  growth,  and  elasto-plastic  coupling  (Maier  and 
Hueckel.  1979). 

2.  Degradation  of  yield  limit,  resulting  from  a  reduction  in  the  areas 
of  cohesive  or  plastic-fractional  connections.  Such  degradation  does 
not  produce  a  reduction  in  the  material  stiffness. 
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The  observed  behavior  of  many  real  materials,  including  concrete  and  rock,  suggests 
that  their  responses  involve  some  combination  of  the  above  two  types  of  degradation. 
The  endochronic  plastic-fracturing  model  described  above  contains  both  types  of 
degradation.  The  softening  function  F  describes  yield  limit  degradation,  while  the 
function  ^  reflects  degradation  in  stiffness. 

With  the  above  discussion  in  mind,  let  us  consider  the  special  cases  of  the 
model  shown  in  Figure  2.  First,  if  the  plasticity  of  the  model  is  suppressed  and  ^  is 
taken  to  depend,  say.  on  the  axial  tensile  strain  e,  then  the  model  reduces  to  the 
classic  elastic-fracturing  model  (see  Dougill.  1976.  for  example),  which  exhibits  only 
stiffness  degradation.  In  this  case,  the  slope  ofthe  unloading/reloeding  path  decreases 
with  increasing  strain,  but  all  unloading  paths  pass  through  the  origin.  The  response 
of  this  special  case  is  shown  in  Figure  2(a).  Secondly,  the  general  model  reduces  to 
the  standard  (non-softening)  endochronic  plasticity  model  when  there  is  no  damage, 
i.e.,  ^  =  F  =  1.  In  this  case,  the  response  of  the  model  is  depicted  in  Figure  2(b). 
Thirdly,  if  one  sets  ^  =  1.  the  general  model  reduces  to  an  endochronic  plastic¬ 
fracturing  model  with  yield  limit  degradation,  whose  response  Is  shown  in  Figure  2(c). 
Note  in  this  case  that  the  model  unloads  and  reloads  with  the  same  slope  as  the 
initial  elastic  slope.  Finally,  when  both  yield  limit  degradation  and  stiffness 
degradation  are  present  (^  ^  1,  F  ^  1),  the  model  has  the  response  features 
illustrated  in  Figure  2(d).  Here,  the  slope  of  the  unloading-reloading  path  decreases 
with  increasing  strain,  and  the  unloading  paths  no  longer  pass  through  the  origin  due 
to  the  development  of  plastic  strains.  These  features  are  representative  of  many 
materials,  including  concrete,  as  we  will  demonstrate  in  the  sequel. 

There  is  also  another  interesting  feature  of  the  model  which  should  be  noted.  If 
the  model  is  subjected  to  monotonic  increasing  tensile  strain,  the  corresponding 
behavior  in  the  transformed  space  (a  -  e)  may  or  may  not  involve  unloading, 
depending  upon  the  relationship  of  the  function  ^  to  a  critical  curve,  which  we  will 
call  This  critical  curve  is  the  curve  in  ^  -  e**  space  which  satisfies  the  condition 
d(^6*’)  =  0.  Given  a  functional  relation  for  say  f  critical  curve  is 

defined  by  the  expression. 


(3-15) 
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Special  cases  of  the  general  endochronic  plastic-fracturing  model  for  simple  tension. 


If  the  curve  lies  above  the  critical  curve  in  the  region  where  d(^e*’)  >  0.  such 

as  curve  A  in  Figure  3.  the  response  in  the  transformed  state  will  involve  only 
loading,  as  shown  in  Figure  4.  This  follows  immediately  from  the  relation  (3-lc).  i.e., 

deP  =  d(^eP]  (3-16) 

On  the  other  hand,  if  the  curve  ^{e^)  falls  below  the  critical  curve  at  some  point  P 
in  the  region  where  d(^e^)  <  0,  such  as  curve  B  in  Figure  3,  unloading  will  take 
place  in  the  transformed  state,  since,  as  shown  in  Figure  4,  in  view  of  Eq.  (3-16), 
de^  will  be  negative.  The  occurence  of  unloading  in  the  transformed  state  does  not, 
however,  affect  the  continuity  or  the  smoothness  of  the  response  in  the  current  state, 
as  our  numer'ical  studies  have  shown.  Thus,  the  model  allows  unloading  (de  <  0)  to 
take  place  in  the  transformed  space  concurrently  with  loading  (de  >  0)  in  the  current 
state. 


In  the  following  section,  an  incremental  numerical  scheme  is  presented  for 
integrating  the  system  of  equations  (3-9)  to  (3-12)  in  conjunction  with  the  evolution 
equation  (3-14), 

3.3  AN  INCREMENTAL  SCHEME  FOR  INTEGRATING  THE 
GOVERNING  EQUATIONS. 

Let  us  now  consider  the  governing  equations  for  the  current  states  which  are 
given  by  Eqs.  (3-9)  to  (3-12)  and  (3-14),  with  E(z)  defined  by  Eq.  (3-8).  If  we 
substitute  Eq.  (3-8)  into  Eq.  (3-9),  we  find  that 


(3-17) 


where 


rZ 


(3-18) 
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Figure  3. 


Figure  4. 


Plot  of  ^  versus  for  two  different  functions  showing  their 

relationship  to  the  critical  curve  a*  (£**). 


Stress  versus  strain  in  the  transformed  space,  showing  the 
corresponding  responses  for  the  two  diffeent  function  depicted 

in  Figure  3. 
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Differentiation  of  Eq.  (3-18)  with  respect  to  z  leads  to  the  expression 

dQ^  +  Q^dz  =  (3-19) 

If  we  now  set: 

n 

(3-20) 

r=l 

n 

r=l 

then,  in  view  of  Eqs.  (3-17)  and  (3-19).  we  can  write: 

da  =  [a  d(^cP]  -  Q  dz]  ♦  ^  d^  (3-22) 

Also.  Eq.  (3-10)  may  be  differentiated  to  give: 

da  =  ^^E^(de  -  de^)  ♦  2^E^[e  -  e^ld^  (3-23) 


Returning  now  to  the  evolution  equation  (3-14).  we  rewrite  this  in  the  form: 
d^  =  H  deP 

where 

H  ^  I"  >  0  ’"<>  '  < 

.0  ,  otherw  i  se 


(3-24) 

(3-25) 


Equation  (3-24)  may  now  be  combined  with  Eqs.  (3-22)  and  (3-23)  to  give  the 
following  expressions: 


da  =  [a^^  ♦  A^e^H  ♦  ^jde^  -  ^Qdz 
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(3-26) 


(3-27) 


<U  =  (^^Ejde  -  2^E^h(€  -  e'>)]deP 


Upon  eliminating  da  from  the  above  equations  and  solving  for  de*^,  we  obtain 
the  expression: 


de^  =  a  de  a.dz 
o  1 


(3-28) 


where 


a  = 
o 


a,  = 


*  E)^  ♦  H(A^eP  -  J] 


(3-29) 


^  (A  ♦  E)^^  ♦  H[A^fP  -  ^ 


From  Eqs.  (3-11),  (3-12)  and  (3-14),  we  can  write 

<)7  •  II  *  "  '"^’‘''’1 


(3-30) 


For  the  case  of  loading,  we  have  de^  >  0,  so  that  Eq.  (3-30)  may  be  written  as 

dz  =  I  ^  1^  ♦  HePl^Oode  ♦  (3-31) 

where  Eq.  (3-28)  has  been  used.  Solving  Eq.  (3-31)  for  dz,  we  obtain 


dz  = 


fl‘r* 

*  HeP 

[-11? 

^  ♦  HeP 

d 

d6 


(de  >  0) 


(3-32) 


Similarly,  for  the  case  of  unloading  (de^  >  0),  we  find 


dz  = 


1  k  -  HdE| 


13  i 


'  •  1!  F  1C  • 


de  (de  <  0) 
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(3-33) 


On  the  basis  of  the  above  developments,  the  following  algorithm  provides  a  numerical 
method  for  incrementally  integrating  the  system  of  equations  (3-9)  to  (3-12)  and 
(3-14)  when  Ae.  a,  e,  F  and  H  are  known  at  time  t; 

1.  Depending  upon  the  sign  of  Ae.  use  either  Eqs.  (3-32)  or  (3-33)  to 
obtain  Az. 

2.  Determine  Ae**  from  Eq.  (3-28). 

3.  Compute  A^  from  Eq.  (3-14), 

4.  Using  Eq.  (3-19),  calculate  AQ^  for  r  =  1,2,  ...  n, 

5.  Evaluate  AQ  from  Eq.  (3-21). 

6.  Determine  La.  using  Eq.  (3-27) 

The  numerical  scheme  described  above  is  first-order  accurate.  It  can  be  made 
second-order  accurate  by  using  Richardson  extrapolation  in  conjunction  with  the  above 
first-order  scheme.  For  further  details,  see  Murakami  and  Read  (1988). 

3.4  APPLICATION  TO  PLAIN  CONCRETE. 

There  are  little  data  presently  available  on  the  complete  load-deformation 
response  of  plain  concrete  to  uniaxial  tension,  and  the  data  that  do  exist  are  often 
conflicting  and  confusing  (compare  the  results  of  Evans  and  Marathe  (1968)  with 
those  of  Petersson  (1981),  for  example).  The  testing  of  brittle,  tension-weak 
materials,  such  as  plain  concrete,  presents  two  major  experimental  problems,  namely, 
(1)  the  tendency  of  a  specimen  to  fail  near  the  grips  where  the  state  of  stress  is  not 
uniaxial  and  (2)  the  inability  of  the  testing  device  to  maintain  stable  post-peak 
response. 

Recently,  Gopalaratnam  and  Shah  (1985)  developed  experimental  techniques  which 
overcome  the  difficulties  noted  above,  and  measured  the  complete  load-deformation 
response  of  plain  concrete  to  simple  tension,  including  the  post-peak  behavior.  Seven 
different  mix  proportions  were  used  to  provide  data  on  the  influence  of  aggregate  size, 
water-cement  ratio,  and  volume  content  of  aggregates.  To  insure  stable  post-peak 
response,  a  small  notch  was  introduced  into  each  specimen  at  its  center  and  on  both 
sides,  and  the  displacement  across  the  notch  was  used  to  control  the  loading. 
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Studies  showed  that  the  average  responses  of  notched  and  unnotched  specimens  were 
identical  and  that  both  notched  and  unnotched  specimens  failed  at  a  single  critical 
section. 

Typical  results  from  the  study  by  Gopalaratnam  and  Shah  (1985)  are  depicted  in 
Figure  5,  where  the  stress-strain  response  of  a  medium  strength  concrete  (f'  = 
6,364  psi)  to  simple  tension  is  shown.  Note  that  the  data  shown  in  this  figure  were 
obtained  from  four  different  types  of  tests.  Three  tests  were  done  under  monotonic 
axial  straining  while  a  fourth  test  was  conducted  under  cyclic  straining  conditions. 
One  test  was  peformed  with  an  unnotched  specimen,  while  the  rest  were  notched. 
The  axial  displacement  for  one  specimen  was  determined  from  optical  measurements  of 
the  crack  width,  while  local  strain  measurements  were  made  on  another  specimen. 
Despite  these  differences  in  testing  procedure,  the  curves  show  a  remarkable 
consistency,  when  one  considers  the  usual  difficulties  of  testing  concrete  in  tension. 

In  the  following  section,  the  endochronic  plastic-fracturing  model  described  in 
Section  3.1  is  applied  to  the  data  of  Gopalaratnam  and  Shah  (1985)  shown  in 
Figure  5.  The  procedures  used  to  determine  the  various  material-dependent 
parameters  and  functions  in  the  model  are  described.  The  response  of  the  resulting 
model  to  monotonic  and  cyclic  tensile  straining  under  uniaxial  loading  is  then 
compared  with  the  corresponding  data. 

3.4.1  Determination  of  Material-Dependent  Parametert  and  Functions. 

The  endochronic  plastic-fracturing  model  for  simple  tension  defined  by  Eqs.  (3-9) 

to  (3-12).  with  an  evolution  equation  for  ^  having  the  form  of  Eq.  (3-14)  will  now  be 

applied  to  Gopalaratnam-Shah  (1985)  data  for  plain  concrete.  As  an  inspection  of 

this  system  of  equations  reveals,  the  material-dependent  functions  and  parameters  that 

must  be  evaluated  from  data  are  E^.  the  Young's  modulus  of  the  undamaged  material; 

E(z).  the  weakly  singular  kernel  function  that  depends  on  the  intrinsic  time,  z;  the 

hardening/softening  function  F;  and  the  evolution  equation  (3-14)  for  The  data 

shown  in  Figure  5  permit  a  direct  determination  of  E  and  4,  but  are  not  sufficient  to 

o  * 

allow  us  to  determine  E(z)  and  F  directly.  Because  of  this,  an  interative  approach 
was  used  to  determine  these  functions  in  the  present  study. 
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From  the  initial  slope  of  the  loading  curve  shown  in  Figure  5,  we  find  that 
=  500  psi.  Figure  6  describes  the  way  in  which  the  elastic  strain  e®.  the  plastic 
strain  e**  and  finally  ^  may  be  obtained  directly  from  the  unloading-reloading  paths 
shown  in  Figure  5.  Note  from  Eq.  (3-10)  that  the  instantaneous  Young's  modulus  E 
at  any  damage  level  is  given  by  E  =  In  this  manner,  the  following  table  was 

constructed  from  the  data  given  in  Figure  5  for  the  five  unloading-reloading  paths. 


Table  I 

Information  on  Current  State  Inferred 
Directly  from  Figure  5. 


a(psi) 

E(psi) 

e'llO-') 

cP(10"*) 

0 

0 

500 

1 

0 

0 

1.4 

295 

357 

0.85 

0.80 

0.60 

2.1 

170 

192 

0.62 

0.90 

1.20 

2.9 

130 

125 

0.50 

1.00 

1.90 

3.7 

95 

81 

0.40 

1.10 

2.60 

4.4 

68 

59 

0.34 

1.20 

3.20 

Using  the  results  given  in  the  above  table,  together  with  Eqs.  (3-1),  the  information 
given  below  in  Table  II  was  derived  for  the  transformed  state. 

Table  II 


Information  on  Transformed  State 
Derived  from  Table  I  Using  Eq.  (3-1). 


e(10-*) 

e(10^ 

a(psi) 

e*(10'^) 

.''(lO-’) 

0 

0 

0 

0 

0 

1.4 

1.18 

349 

0.68 

0.51 

2.1 

1.30 

27 

0.56 

0.77 

2.9 

1.45 

260 

0.50 

1.00 

3.7 

1.48 

237 

0.44 

1.04 

4.4 

1.51 

198 

0.41 

1.13 

The  function  f  describes  the  degradation  of  material  stiffness,  caused  by 
microcracking,  void  growth  and  elastic-plastic  coupling  (Maier  and  Hueckel,  1979).  As 
noted  earlier,  we  found  that  ^  is  most  naturally  described  in  terms  of  the  plastic 
strain  e**.  Guided  by  the  data  given  in  Table  I,  ^  was  taken  in  the  form: 
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Average  Net  Axial  Stress,  o(psi) 


eoo 


Axial  Strain  (tistr) 


Figure  5.  Response  of  plain  concrete  to  uniaxial  tension,  as  measured  in  four 
different  types  of  tests  (from  Gopolaratnam  and  Shah,  1985). 


and  ^  are  defined  at  an  arbitrary  unloading  point  P. 
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^  =  1  -  [1  -  exp 


(3  34) 


Inhere 


k  =  9.33  X  10" 


m  =  1.21 


K  =  0-20 


(3-35) 


Figure  7  depicts  this  relationship  and  demonstrates  its  ability  to  describe  the  data. 

With  ^  defined,  the  evolution  equation  for  ^  is  then  simply  obtained  by 
differentiating  Eq.  (3-34).  which  leads  to  the  result 


=  H(cP]deP 


where 


(3-36) 


(3-37) 


(3-38) 


Equation  (3-36).  with  H(e^)  defined  by  Eq.  (3-37)  together  with  the  values  of  the 
parameters  given  in  Eq.  (3-35).  is  the  evolution  equation  for  f  used  in  the  present 
study. 

Turning  now  to  the  kernel  function  E(z)  and  the  softening  function  F,  we  note 
that  they  are  both  associated  with  the  plasticity  of  the  model  and.  as  such,  are  most 
naturally  defined  with  respect  to  the  response  in  the  transformed  state,  particularly, 
the  response  in  the  a  -  plane.  We  note  that  the  response  of  the  model  in  the 
ff  -  c**  plane  is  governed  by  Eqs.  (3-2).  (3-4)  and  (3-5).  Assuming  that  F  depends 
upon  6^.  these  equations  may  be  combined  to  give  the  expression: 


I 


•)  f[£P(z-)1<Iz' 


(3-39) 


I 


Dependence  of  ^  on  plastic  Figure  B.  Stress  versus  plastic  strain  in  the 

strain  for  monotonic  extension  transformed  state  for  monotonic 

under  simple  tension.  extension  under  simple  tension,  showing  yield 

limit  degradation. 


where 


fl  Ide^l 
fG") 


(3-40) 


anci  the  unknown  functions  E(z)  and  F{£*’)  are  to  be  determined.  As  Eq.  (3-39) 
reveals,  these  functions  appear  in  product  form,  which  adds  complexity  to  their 
determination. 


As  reference  to  Table  II  will  show,  the  data  defining  the  behavior  in  the  a  -  ^ 
plane  are  very  sparse  (5  points)  and  there  are  no  data  points  for  small  values  of 
where  information  is  required  if  E(z)  is  to  be  evaluated  directly.  As  a  result,  it  was 
necessary  to  determine  E(z)  and  F(6*’)  indirectly,  using  an  iterative  approach. 
Proceeding  in  this  manner,  the  following  expressions  for  E(z)  and  F(e*’)  were  found  to 
describe  the  behavior  of  the  concrete  quite  well: 


E(2) 


(3-41) 


where 


Aj  =  1.62  X  10^  psl 
A^  =  2.69  X  10*  psi 


Oj  =8.82  X  10^ 
Oj  =  3.99  X  10® 


(3-42) 


and 


F(eP}  =  F.  >  (1  -  F.j  exp[-  aleP)"] 


(3-43) 


where 

fl  =  6.71  X  10^ 

n  =0.866  (3-44) 

F„  =  0.20 

«P 

Based  on  the  forms  for  E(z)  and  F(6P)  given  above,  the  response  of  the  model  in  the 
transformed  {a  -  e^)  space  to  monotonic  uniaxial  loading  is  given  in  Figure  8.  The 
corresponding  data  points  from  Table  II  are  also  shown  in  this  figure  for  comparison. 
As  the  figure  shows,  the  model  describes  the  data  quite  well  in  this  space. 


B-35 


The  specification  of  the  paranrieters  and  functions  in  the  model  is  now  complete. 


3.4.2  Predicted  and  Measured  Responses. 


Using  the  parameters  and  functions  given  above,  response  of  the  endochronic 
plastic-fracturing  model  defined  by  Eqs.  (3-9)  to  (3-12)  was  examined  for  both 
monotonic  and  cyclic  tensile  strain  paths  over  the  range  of  tensile  strains  studied  by 
Gopalaratnam  and  Shah  (1985).  The  governing  equations  of  the  model  were 
integrated  using  the  numerical  scheme  described  in  Section  3.3.  The  results  obtained 
in  this  manner  are  depicted  in  Figure  9.  where  both  the  model  predictions  and  the 
corresponding  data  are  shown  for  comparative  purposes.  As  the  figure  reveals,  the 
model  describes  the  data  quite  well.  The  only  noticeable  difference  between  the 
predicted  and  observed  responses  is  that  the  data  show  significant  hysteresis  during 
the  unloading-reloading  processes  while  the  model  does  not.  Overall,  however,  the 
agreement  is  considered  very  encouraging. 

Finally,  we  note  that  during  the  course  of  calculating  the  response  to  monotonic 
straining  shown  in  Figure  9,  the  dependence  of  ^  on  both  the  total  strain  e  and  the 
elastic  strain  c*  was  monitored.  The  resulting  curves  are  shown  in  Figure  10,  where 
it  is  seen  that  ^  decreases  in  a  monotonic  manner  with  e.  but  is  a  double-valued 
function  of  e*.  This  clearly  shows  why  attempts  to  define  ^  as  a  simple  monotonic 
decreasing  function  of  €*  proved  futile. 
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Average  Net  Axial  Stress,  o(psi) 


0  100  200  300  400 


Axial  Strain  (/latr) 

(a)  Response  measured  by  Gopalaratnam  and  Shah  (1985). 


100  200  300  4 

Axial  Strain  (^str) 

(b)  Response  predicted  by  present  model. 


Figure  9.  Response  of  plain  concrete  to  simple  tension. 
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Figure  10.  Dependence  of  ^  on  £  and  £*,  as  predicted  by  the  present  constitutive  model  for  plaii 
concrete. 


Section  4 


CONCLUSIONS 

This  study  represents  the  first  attempt  to  explore  the  application  of  the  new 
endochronic  plastic-fracturing  theory  to  plain  concrete.  For  this  purpose,  attention  was 
focussed  on  the  behavior  of  the  model  under  uniaxial  loading  conditions.  On  the 
basis  of  the  results  presented  herein,  the  following  conclusions  are  drawn. 

1.  To  the  extent  that  the  model  has  been  explored  in  the  present 
study,  i.e..  under  conditions  of  simple  tension,  the  model  appears 
to  provide  an  excellent  description  of  the  observed  behavior  of  plain 
concrete  during  progressive  failure  under  tensile  straining  conditions. 

2.  In  addition  to  elasto-plastic  effects,  the  model  can  describe  both 
stiffness  degradation,  through  and  yield  limit  degradation 
through  F. 

3.  The  endochronic  plastic-fracturing  model  contains,  as  special  cases, 
several  models  that  are  noteworthy.  This  includes  the  classic 
elastic-fracturing  model,  the  standard  plasticity  model  and  a  plastic¬ 
fracturing  model  with  yield  degradation.  The  general  model  includes 
both  stiffness  degradation  and  yield  limit  degradation. 

4.  The  key  to  the  success  of  the  model  lies  in  having  the  appropriate 
evolution  equation  for  ^  for  the  material  of  interest.  Ideally,  it 
would  be  desirable  to  have  the  evolution  equation  reflect  the 
underlying  micromechanical  fracture  mechanisms,  but  the  technology 
has  not  advanced  to  that  point.  A  research  effort  to  develop  such 
an  evolution  equation  for  ^  is  currently  underway  (Valanis,  1988), 
but  the  results  from  this  study  are  not  expected  to  be  available  for 
sometimes. 

5.  In  the  present  study,  an  evolution  equation  for  ^  suggested  by  the 
data  was  used.  We  found  that  this  equation  was  most  naturally 
expressed  in  terms  of  the  plastic  strain  e**.  The  dependence  of  ^ 
on  the  elastic  strain  was  found  to  be  double-valued. 


B-39 


6. 


A  simple  and  efficient  numerical  scheme  was  developed  for 
integrating  the  governing  equations  of  the  model  for  simple  tension, 
and  this  suggest  an  approach  for  numerically  dealing  with  the  more 
general  three-dimensional  version  of  the  model. 


7.  Procedures  were  developed  for  determining  the  materials-dependent 
parameters  and  functions  in  the  model  from  standard  data. 
However,  because  of  the  sparseness  of  the  available  data,  several  of 
the  functions  had  to  be  determined  indirectly. 


8. 


On  the  basis  of  the  very  encouraging  results  presented  herein, 
further  studies  should  be  undertaken  to  explore  and'  validate  the 
model  under  more  general  loading  conditions,  such  as  uniaxial 
compression  and  triaxial  compression.  These  cases  involve  cracking 
patterns  that  are  considerably  more  complex  than  occurs  under 
simple  tension  and  thus  should  provide  a  stringent  test  of  the 
theory. 
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